NAVAL  POSTGRADUATE  SCHOOL 

Monterey,  California 


AD-A208  536 


DEVELOPMENT  OF  A  SYSTEM  MODEL  AND 
LEAST  MEAN  SQUARE  (LMS) 

FILTER  FOR  THE  NAVAL  POSTGR.4DUATE 
SCHOOL  (NPS) 

INFRARED  SEARCH  AND  TARGET  DESIGNATION 
SYSTEM  (IRSTD) 

by 

Michael  Louis  Gribaudo 
March  19S9 


Thesis  Advisor  Alfred  W.  Cooper 


Approved  for  public  release;  distribution  is  unlimited. 


89605  121 


Unclassified 


secLiriu  c!assif;cat!an  of  this  pa;;e 


la  Report  Security  Ciassificat'oi-  Unclassified 


2a  Security  Clasiification  Authi'rity 


dine  Schedule 


non  Rcpo't  \urr,beris) 


ha  .Name  of  Performing  Organization 
Naval  Postgraduate  School 


cc  Address  icin',  siuie.  anj  ZIP  cede' 

Monterev.  CA  93943-5000 


Sa  Name  of  Funding  Sponsoring  Organization 


REPORT  DOCUMENTATION  PAGE 


lb  Restriciise  .Markings 


6  b  Office  Symbol 
( if  apj'iicable}  6 1 


8b  Office  Symbol 
{if  applicable) 


3  Distribution  Availability  of  Report 

Approved  for  public  release:  distribution  is  unlimited. 


5  Monitoring  Organization  Report  Number(s) 


7a  .Name  of  Monitoring  Organization 

Naval  Postgraduate  School 


7b  Address  (city,  state,  and  ZIP  code) 

Monterev.  CA  93943-5000 


9  Procurement  Instrument  Identification  .Number 


'.  and  ZIP  cede,' 

\  10  Source  of  Funding  Numbers 

1  Program  Element  No  I 

1  Project  No  I 

[Task  No  1 

1 1  Title  {in  iude  securir:  cia.';.\iif  alien  <  DEVELOPMENT  OF  A  SYSTEM  MODEL  AND  LEAST  MEAN  SQUARE  (LMS) 
FILTER  FOR  THE  N.AV AL  POSTGRADUATE  SCHOOL  (NPS)  INFRARED  SEARCH  AND  TARGE!  DESIG- 
N.ATION  (IRSTD)  SYSTEM 


12  Personal  .Author' si  Michael  Louis  Gribaudo 


13a  Type  of  Report 

Ma?ter  s  Thesis 


13b  1  ime  Covered 
From  To 


Id  Date  of  Report  (year,  month,  dayi 

March  1989 


1 5  Paae  Count 
108 


16  Supplementary  Notation  !  he  views  expressed  in  this  thesis  are  those  of  the  author  and  do  not  reflect  the  official  policy  or  po¬ 
sition  of  the  Depiulmcnt  of  Defense  or  the  U.S.  Government. 


C  .,'ues 


Group 


l.s  Subject  Terms  (continue  on  reverse  if  necessary  and  identity  by  ble.k  nutnbei  i 

IRS  1  D.  LMS  ftlter,  S  N  ratio.  Signal  Processing 


19  .-Abstract  iccr.p.nue,  on  rcstrse  v  r.f'  cssary  and  identify  by  block  number) 

^  .A  system  model  and  a  least  mean  square  (LNIS)  filter  for  the  Naval  Postgraduate  School  (NPS)  Infrared  Search  and  Target 
Designation  (IRS  I  D)  system  were  developed. 

!  he  system  model  was  developed  and  run  on  the  NPS  IBM  3033  4381  mainframe  computer  network.  The  model  sim¬ 
ulated  the  eflects  of  the  optics  and  electronic  processing  equipment  of  the  IRSTD  system,  and  produced  output  data  repre¬ 
sentative  of  the  detector  outputs  of  the  system. 

The  outputs  of  the  IRS  I  D  model  were  used  to  develop  a  digital  filter  based  on  the  principle  of  least  mean  square  opti¬ 
mization  between  an  actual  IRSTD  detector  output  and  a  power  series  expansion  representing  a  detector  output  containing 
both  background  clutter  and  a  model  target  signal.  It  was  determined  that  the  raised  cosine  function  ser\ed  as  the  best  model 
for  IRSTD  point  and  near-point  targets,  (0  1  mrad  by  0.1  mrad  to  1.5  mrad  by  1.5  mrad),  and  a  set  of  trial  L.MS  filters  were 
generated  based  on  this  model. 

After  filtering  both  simulated  and  real  data,  consisting  of  simulated  and  real  target  signals  embedded  in  simulated  and  real 
backgrounds,  it  was  determined  that  an  LMS  filter  generated  from  a  raised  cosine  with  a  half-amplitude  width  of  0.9  mrad 
was  optimal  for  point  and  near-point  targets.  The  signal-to-noise  ratios  of  all  target  and  background  combinations  increased 
by  a  factor  of  appro.\imatel\  30  for  the  simulated  backgrounds,  and  approximately  six  for  the  real  backgrounds,  upon  filtering 
the  detector  outputs  with  the  optimal  L.MS  filter,  f 

It  is  believed  that  this  filter  should  be  incorporated  into  the  NPS  IRSTD  system  as  an  initial  signal  processing  filter,  and 
that  the  filtered  outputs  are  appropriate  for  use  as  inputs  to  target  detection  and  acquisition  routines. 


20  Distribution  .Availability  of  Abstract 
S  unclassified  unlimited  D  same  as  report 


22a  .Name  oi  Responsible  Individ'jal 

.Alfred  W.  Cooper 


DD  FORM  1473.84  MAR 


21  Abstract  Security  ClassiFication 
□  DTic  users  Unclassified 


22b  Telephone  (include  Area  coder 

(408)  646-2452 


8.3  APR  edition  may  be  used  until  exhausted 
All  other  editions  are  obsolete 


22c  Office  Svmbol 

6ICr 


security  classifcaiion  of  this  page 

I  nclassified 


1 


Approved  for  public  release;  distribution  is  unlimited. 


Development  of  a  System  Model  and  Least  Mean  Square  (LMS) 
Filter  for  the  Naval  Postgraduate  School  (NPS) 

Infrared  Search  and  Target  Designation  (IRSTD)  System 

by 

Michael  Louis  Gribaudo 
Lieutenant,  United  States  \a\7 
B.S.,  University  of  Idaho,  1983 

Submitted  in  partial  fulfillment  of  the 
requirements  for  the  degree  of 

MASTER  OF  SCIENCE  IN  PHYSICS 

from  the 

NAVAL  POSTGRADUAl  E  SCHOOL 
March  1989 


ABSTRACT 


A  system  model  and  a  least  mean  square  (LMS)  filter  for  the  Naval  Postgraduate 
School  (NPS)  Infrared  Search  and  Target  Designation  (IRSTD)  system  were  developed. 

The  system  model  was  developed  and  run  on  the  NTS  IBM  3033/4381  mainframe 
computer  network.  The  model  simulated  the  effects  of  the  optics  and  electronic  proc¬ 
essing  equipment  of  the  IRSTD  system,  and  produced  output  data  representative  of  the 
detector  outputs  of  the  system. 

The  outputs  of  the  IRSTD  model  were  used  to  develop  a  digital  filter  based  on  the 
principle  of  least  mean  square  optimization  between  an  actual  IRSTD  detector  output 
and  a  power  series  expansion  representing  a  detector  output  containing  both  back¬ 
ground  clutter  and  a  model  target  signal.  It  was  determined  that  the  raised  cosine 
function  served  as  the  best  model  for  IRSTD  point  and  near-point  targets,  (0.1  mrad  by 
0.1  mrad  to  1.5  mrad  by  1.5  mrad),  and  a  set  of  trial  LMS  filters  were  generated  based 
on  this  model. 

.After  filtering  both  simulated  and  real  data,  consisting  of  simulated  and  real  target 
signals  embedded  in  simulated  and  real  backgrounds,  it  was  determined  that  an  LMS 
filter  generated  from  a  raised  cosine  with  a  half-amplitude  width  of  0.9  mrad  was  optimal 
for  point  and  near-point  targets.  The  signal-to-noise  ratios  of  all  target  and  background 
combinations  increased  by  a  factor  of  approximately  30  for  the  simulated  backgrounds, 
and  approximately  six  for  the  real  backgrounds,  upon  filtering  the  detector  outputs  with 
the  optimal  L.MS  filter. 

It  is  believed  that  this  filter  should  be  incorporated  into  the  NPS  IRSTD  system  as 
an  initial  signal  processing  filter,  and  that  the  filtered  outputs  are  appropriate  for  use  as 
inputs  to  target  detection  and  acquisition  routines. 
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1.  INTRODUCTION 


A.  GENERAL 

Several  years  after  completion  of  at-sea  trials  in  L'SS  Kinkaid,  the  Advanced  De¬ 
velopment  Model  (ADM)  of  the  SPAR  Aerospace  AN  SAR-8  Infrared  Search  and 
Target  Designation  (IRSTD)  system  was  transferred  to  the  Naval  Academic  Center  for 
Infrared  Technology  (NACIT)  at  the  Naval  Postgraduate  School  (NPS)  in  Monterey. 
California.  NACI  T  received  this  system  with  the  intention  of  using  it  in  support  of  re¬ 
search  in  infrared  physics,  infrared  technology,  and  infrared  signal  processing.  Specific 
objectives  of  this  research  include: 

•  Collection  and  study  of  unfiltered  infrared  background  scenes. 

•  Development  of  signal  processing  techniques  for  the  rejection  of  background  clut¬ 
ter. 

•  Dcw'lopment  of  target  detection,  acquisition,  tracking,  and  classification  tech¬ 
niques. 

Since  receiving  the  IRSI  D  system.  NACIT  personnel  have  accomplished  a  great 
deal  of  work  involving  the  system  installation,  testing,  and  callibration.  While  installing 
the  system  in  Spanagal  Hall  at  NPS.  a  number  of  modifications  were  required  to  be 
made.  Some  of  the  more  important  modifications  (Ref.  1)  include; 

•  Removal  of  the  Stable  Platform  and  Stable  Platform  Control  Console. 

•  Bx  pass  of  the  Background  N'ormali/cr. 

•  Replacement  of  the  Data  Processing  Unit  with  a  MASSCOMP  computer  system 
of  higher  capacify. 

•  Replacement  ol'  the  cryogenic  engine  and  insulating  vacuum  assembly,  leading  to 
the  removal  of  the  inner  germanium  window  of  the  system's  optical  subassembly 
(OSA). 

Of  these  modifications,  the  ciy  ogenic  engine  replacement  had  the  most  impact  on 
the  operation  of  the  system.  Removal  of  the  germanium  window  changed  the  optical 
path  of  the  system,  thereby  increasing  the  system  s  optical  spot  size.  To  correct  for  this 
eflect,  it  was  found  [Refs.  1.2]  that  a  focal  shift  of  the  detector  array  by  approximately 
0.133  inches  towards  the  reflecting  mirror  restored  the  modified  system  to  a  level  of 
performance  comparable  to  that  of  the  original  system. 
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Upon  restoring  the  system  to  an  acceptable  level  of  performance,  the  relative 
responsivities  of  the  individual  detector  elements  were  calibrated  with  their  associated 
electronics.  It  was  Ibund  [Ref  2]  that  seven  lead  array  detectors  and  eight  lag  array 
detectors  were  non-functional,  and  that  three  additional  lag  array  detectors  were  per¬ 
forming  at  levels  below  standards.  It  was  also  found  that  the  relative  response  ampli¬ 
tudes  of  the  lag  array  detectors  were  a  factor  of  1.85  higher  than  the  lead  array  detectors 
for  the  same  calibration  source.  This  dinerence  was  attributed  to  the  use  of  din'erent 
bandwidth  filters  for  the  two  detector  arrays. 

As  of  January  1989.  The  system  is  in  place  in  Spanagel  Hall  at  NPS  .Monterey,  and 
is  ready  to  take  and  store  data  in  the  MASSCOMP  computer.  With  these  capabilities 
available,  the  N,-\CIT  research  group  can  begin  collecting  infrared  background  data,  and 
can  also  begin  developing  clutter  rejection  filters. 

B.  PI  RPOSE  AND  ORGANIZATION  OF  THESIS 

The  purpose  of  this  thesis  is  to  develop  a  digital  filter  capable  of  extracting  point 
and  near-point  targets  from  typical  infrared  background  clutter  scenes.  The  filter  se¬ 
lected  for  this  task  is  constructed  using  the  principle  of  Icast-mean-square  optimization, 
and  is  tlierefore  known  as  an  LMS  filter  [Ref  3]. 

The  text  of  this  thesis  describes  the  derivation  of  the  LMS  filter,  and  presents  the 
results  of  its  use  on  both  simulated  and  real  data.  The  thesis  is  organized  as  follows: 

•  Chapter  II  presents  a  brief  description  of  the  NPS  IRSTD  system,  with  particular 
emphasis  on  the  ss  stem's  specifications. 

•  Chapter  Ill  discusses  tlic  nature  of  infrared  backgrounds,  as  compared  to  targets 
o!  interest  to  IRSTD  signal  processing,  and  also  presents  an  inlVared  background 
simulation  model. 

•  Chapter  IV  describes  the  dcvclopemcnt  of  an  IRSTD  simulation  model,  con¬ 
structed  to  generate  simulated  IRSTD  output  signals  for  initial  analysis  of  the  LMS 
filter. 

•  Chapter  V  describes  the  developement  of  an  LMS  filter  for  extracting  targets  of 
interest  from  infrared  background  scenes. 

•  Chapter  VI  presents  the  results  of  filtering  both  simulated  and  real  infrared  IRSTD 
data  with  the  LMS  filter. 

•  Chapter  VII  discusses  the  results  of  chapter  VI  and  presents  recommendations  for 
further  studies  in  this  area. 

•  .Appendix  A  presents  plots  of  filtered  outputs  illustrating  the  increase  in  target  and 
noise  amplitude  with  increasing  filter  width. 

•  Appendix  B  presents  plots  of  filtered  outputs  of  all  target  and  background  combi¬ 
nations  studied  in  this  thesis,  filtered  with  their  proposed  optimal  L.MS  filters. 


Appendix  C  lists  the  background  simulation  model  FORTRAN  code. 
Appendix  D  lists  the  IRSTD  simulation  model  FORTRAN  code. 
Appendix  E  lists  the  LMS  Filter  FORTRAN  code. 


II.  THE  N PS  IRSTD  SYSTEM 


A.  SYSTEM  DESCRIPTION  AND  SET-UP 

The  NFS  IRSTD  system  is  composed  of  four  major  subsystems.  Since  only  a  brief 
description  of  these  subsystems  will  be  presented  in  this  thesis,  the  reader  is  directed  to 
the  system  technical  manual  [Ref.  4],  for  a  detailed  description  of  the  system. 

A  brief  description  of  the  function  and  set-up  of  the  major  subsystems  and  various 
peripheral  components  of  the  NFS  IRSTD  system  follows; 

•  The  Scanner  Assembly,  located  on  the  roof  of  Spanagel  Hall  at  NFS  .Monterey, 
houses  the  Optical  Subassembly  (OS.A).  The  OSA  consists  of  the  system  optics, 
the  lead  and  lag  detector  arrays,  and  the  cryogenics.  As  a  unit,  the  Scanner  As¬ 
sembly  physically  rotates  the  OS.A  in  complete  360  degree  scans.  The  OSA,  in  turn, 
collects  and  responds  to  infrared  radiation  incident  on  the  OSA  aperture.  Elec¬ 
tronic  signals  from  the  detector  arrays  are  then  transmitted  to  the  Data  Condi¬ 
tioner  Unit  via  the  Buffer  Fower  Unit. 

•  The  Scanner  Assembly  Control  Console  (SACC),  located  on  the  seventh  floor  of 
Spanagel  Hall,  supplies  power  to,  and  controls,  the  Scanner  .Assembly  and  its  as¬ 
sociated  electronic  components.  The  SAAC  controls  parameters  such  as  the 
Scanner  .Assembly  scan  rate. 

•  The  Data  Conditioner  Unit  (DCU),  also  located  on  the  seventh  floor  of  Spanagel 
Hall,  receives  signals  from  the  OS.A  via  the  Buffer  Fower  Unit.  It  is  in  the  DCU 
that  analog-io-digital  conversion,  multiplexing,  and  pre-processing  of  the  detector 
outputs  occurs. 

•  The  .Vl.ASSCO.VIF  computer,  located  on  the  second  floor  of  Spanagel  Hall.  recei%es 
digital  detector  output  signals  from  the  DCU  for  processing,  display,  and  storage. 

B.  SYSTEM  SPECIFICATIONS 

In  order  to  proceed  effectively  with  the  development  of  an  IRSTD  model  and  the 
LMS  filter,  it  is  essential  to  know  some  of  the  IRSTD  system's  technical  specifications. 
.A  short  listing  of  the  general  specifications  is  presented  below  [Refs.  1,2,  and  4]: 

•  Number  and  type  of  detector  elements  -  90  indium  antimonide  (InSb)  detector  el¬ 
ements  per  array  (2  arrays). 

•  Size  of  detector  elements  -  2.0  by  0.3  milliradians  (0.3  milliradians  along  the  scan 
dimension). 

•  OSA  scan  rate  -  variable  (designed  to  operate  at  one-half  revolution  per  second). 

•  DCU  sampling  rate  -  10,000  8-bit  samples  per  second;  the  8-bit  analog-to-digital 
converter  will  soon  be  replaced  by  a  12-bit  converter.  (This  sampling  rate  corre¬ 
sponds  to  0.10472  mrad  per  sample  at  one-half  revolution  per  second.) 

•  Detector  dwell  time  -  approximately  3  samples  per  detector  dwell. 


4 


•  optical  spot  size  -  approximately  0.7  mrad  in  diameter,  based  on  the  diameter  of 
the  blur  spot  of  a  point  target  at  the  detector  array. 

•  System  bandwidth  filters  -  3  to  5  micrometers. 

•  Optics  -FI  Schmidt  catadioptric  telescope. 

•  OSA  aperture  -  10  inch  diameter  with  germanium  aspheric  corrector. 

•  Focusing  mirror  -  16.625  inch  diameter  aluminum  spherical  mirror. 

•  Coolant  -  liquid  nitrogen,  resulting  in  a  detector  operating  temperature  of  83K. 
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III.  INFRARED  BACKGROUNDS 


A.  THE  NATURE  OF  INFRARED  BACKGROUNDS 

Typical  infrared  backgrounds  consist  of  scattered,  reflected,  or  emitted  infrared  ra¬ 
diation  from  objects  other  than  those  sought  by  the  system.  These  non-target  objects 
include,  but  are  not  restricted  to; 

•  the  sun 

•  clouds 

•  precipitation 

•  sea  spray 

•  land 

•  birds 

The  targets  of  interest,  for  the  purpose  of  this  thesis,  are  restricted  to  aircraft  and 
missiles.  Since  the  spatial  extent  of  these  targets  is  typically  ver\'  small,  they  can  be 
modeled  as  point,  or  near-point,  targcts.l  It  is  therefore  only  necessarj’  to  consider  how 
the  spatial  extent  of  these  targets  difiers  from  the  spatial  extent  of  the  background  fea¬ 
tures. 

The  statistical  properties  of  infrared  backgrounds  are  available  from  previous  studies 
[Ref  6].  These  studies  show  that  a  statistical  model  of  background  radiance  can  be  re¬ 
garded  as  a  random  set  of  two-dimensional  pulses  whose  amplitudes  can  be  described 
by  either  Gaussian  or  Poisson  distributions,  depending  on  the  wavelength  and  type  of 
infrared  radiation,  and  whose  pulse  widths  can  be  described  by  Poisson  distributions. 

In  the  3-5^m  region,  the  radiance  received  by  the  system  is  due  to  a  combination 
of  scattered  sunlight  and  thermal  radiation  emitted  from  the  background.  Since  the 
amplitude  of  the  radiance  of  scattered  sunlight  obeys  Poisson  statistics,  and  the  ampli¬ 
tude  of  the  radiance  of  thermal  radiation  obeys  Gaussian  statistics,  the  amplitude  of  the 

1  For  example,  a  typical  jet  aircraft  engine  has  an  exhaust  nozzle  with  an  area  approximately 
3600  cm^  [Ref.  5|.  which  corresponds  to  a  diameter  of  approximately  68  cm.  At  a  range  of  10 
kilometers  (km)  this  source  will  subtend  an  angle  of  0.068  milliradians  (mrad),  at  5  km  it  subtends 
0.136  mrad,  at  1  km  it  subtends  0.68  mrad,  and  at  500  meters  it  subtends  1.36  mrad.  Therefore, 
typical  air  targets  can  be  modeled  as  targets  which  subtend  angles  ranging  from  0.1  to  1.5  mrad, 
where  due  to  the  sampling  rate  of  the  DCU  a  0.1  mrad  target  is  considered  to  be  a  point  target. 


6 


radiance  of  the  background  in  the  3-5Mm  region  should  be  describable  by  some  combi¬ 
nation  of  Gaussian  and  Poisson  distributions  [Ref.  6]. 

Figure  1  is  a  plot  of  the  amplitude  distribution  of  the  radiance  of  a  sky  background 
in  the  3-5.um  region,  as  determined  by  experimental  observation  [Ref.  6].  It  is  interesting 
to  note  that  the  influence  of  the  scattered  sunlight  only  slightly  distorts  the  amplitude 
distribution  from  a  true  Gaussian  form.  Additionally,  it  should  be  noted  that  the  mean 
value  of  the  background  radiance  is  greater  than  zero.  However,  the  mean  radiance  can 
be  changed  to  zero,  or  any  other  value,  by  adding  or  subtracting  a  constant  from  the 
individual  radiances  measured  over  the  entire  background.  This  is  often  desirable,  when 
plotting  background  radiance  over  a  region  of  space,  because  it  can  be  used  to  change 
the  offset  of  the  plot  while  leaving  all  other  statistical  properties  unaltered. 


CD  ^ 


rmplitude: 


Figure  1.  Amplitude  Distribution  of  Background  Radiance:  Solid  line,  background 
radiance  in  the  3-5(um  region.  Dotted  line,  Gaussian  distribution  over- 
laied  for  comparison.  Radiance  units  normalized  to  a  peak  amplitude 
of  ±  0.5.  [Ref.  6] 
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In  the  following  section,  the  development  of  a  background  radiance  distribution  in 
the  spatial  domain  will  be  presented.  In  order  to  simplify  calculations,  the  amplitude 
distribution  of  the  background  radiance  will  be  modeled  as  a  Gaussian  function,  having 
a  mean  value  of  zero  and  a  standard  deviation  ranging  from  0.1  to  0.15.  The  back¬ 
grounds  generated  from  this  model  w'ill  then  be  statistically  analyzed  to  ensure  that  their 
radiance  distributions  are  representative  of  the  distribution  shown  in  Figure  1  on  page 
7;  then  they  will  be  used  to  simulate  IR  backgrounds  for  the  purpose  of  evaluating  the 
effectiveness  of  the  proposed  LMS  filter. 

B.  THE  BACKGROUND  MODEL 

A  background  model  for  the  3-5A<m  region  is  developed  using  the  statistical  proper¬ 
ties  outlined  in  the  previous  section.  The  amplitude  distribution  of  the  background  ra¬ 
diance  can  be  described  by  the  probability  density  function 

P(B:.)  - - =L=^  e"  (3-1) 

where  B,,  is  an  element  of  the  background  radiance  matrix  B  which  describes  the  radiance 
at  given  points  in  the  background  distribution,  B  is  the  average  radiance  over  the  entire 
background,  and  cr*  is  the  variance  of  the  radiance  over  the  entire  background. 

.Adjacent  points  in  the  background  distribution  are  spatially  correlated,  with  a  cor¬ 
relation  function  of 

L(r)  =  (3  -  2) 

where  a"'  is  the  correlation  length,  and  r  is  the  distance  between  any  two  adjacent  points. 

Following  the  model  developed  by  Ben- Yosef  eial.  [Ref.  7  J,  the  background  radi¬ 
ance  distribution  can  be  generated  directly  in  the  space  domain.  This  is  accomplished 
by  determining  the  conditional  probability  density  function,  which  describes  the  proba¬ 
bility  of  finding  if  the  radiances  of  adjacent  points  arc  known. 
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As  a  first  approximation,  B  can  be  taken  as  zero,  and  the  conditional  probability 
density  function  for  determining  the  radiance  B,„  from  the  radiance  B,*  or  B^,,  where 
k=J  —  I  and  /  =  /  —  1,  can  be  written  as 


P(BylB,,)  = 


J _  -(By-ZBj 


y.2nafi^ 


(3-3) 


where 


Z  =  e‘ 


(3-4) 


and 


(3-5) 


Likewise,  the  conditional  probability  density  function  for  determining  the  radiance 
B,.  from  the  radiances  B,*  and  B,,  is 


27r<7/>2 


ZB,, 


(!+Z^^^)  d+Z'-'^) 


ZB;  .  ,  , 

- ^  '  2a 


(3-6) 


where 


/?2  = 


(1  -2Z^) 


(5-7) 


V  (1+Z^") 

From  Equation  (3-3),  it  follows  that  given  B,*,  B,,  is  normally  distributed  with  a  mean 


of 


<B^.>=ZB;fe  (3-8) 

and  a  standard  deviation  of  aP^.  Additionally,  Equation  (3-6)  shows  that  given  B,*  and 
B,„  B„  is  normally  distributed  with  a  mean  of 

<Bp  - - ^^(B,*-f  B^.)  (3-9) 

I  -I- 


and  a  standard  deviation  of  o/Jj. 
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The  background  scene  can  therefore  be  generated  using  any  normally  distributed 
random  number  generation  scheme.  For  this  thesis,  the  International  Mathematical  and 
Statistical  Libraries  (IMSL)  subroutine  GGXML  (Ref.  8]  was  used  to  generate  a  set  of 
normally  distributed  random  numbers,  with  a  mean  of  zero  and  a  standard  deviation  of 
one.  The  random  numbers  were  then  placed  in  all  but  the  first  element  of  a  4000  by  20 
element  matrix.  The  elements  of  this  matrix  can  be  represented  as  Rjj,  with  the  index  i 
ranging  from  one  to  4000  and  the  index  j  ranging  from  one  to  20, 

The  elements  of  the  background  scene  (B,^)  are  then  generated  in  a  matrix  of  similar 
dimensions  by  first  setting  the  value  of  the  origin  equal  to  zero;  the  average  radiance  B. 
The  remaining  elements  of  the  first  row  and  first  column  are  then  determined  using  the 
results  of  Equation  (3-3),  and  can  be  represented  as 

B,-i  =  <T/?i  Ra  -f  <  B„  >  (2  <  /  <  4000)  (3  -  10) 

for  the  first  row.  and 

B,j  =  oP,Ry+<By>  {2<j<20)  (3-11) 

for  the  first  coluinn.  From  Equation  (3-8).  the  mean  values  <B,i>  and  <B,^>  repre¬ 
sent 

<B„>=ZB,,  (3-12) 

and 

<B,y>=ZB„  (3-13) 

In  a  similar  manner,  the  remaining  elements  of  the  background  radiance  distribution 
can  be  determined  from  the  results  of  Equation  (3-6)  and  are  represented  as 

B,y  =  oP2^ij  +  <  By  >  (3-14) 

where  2  <  /  ^  4000  and  2^j^  20.  From  Equation  (3-9),  the  mean  value  <  B,>  is  re¬ 
presented  by 

<B„> - +  (3-15) 

(1  -l-Z^^) 
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It  follows  from  this  background  generation  scheme,  that  the  resultant  background 
radiance  distribution  for  a  background  with  B  =  0  is  dependent  only  on  the  two  variables 
o*  and  (x  ',  and  the  random  number  generator  seed  which  produces  a  different  set  of  Ry 
for  each  seed  used.  Additionally,  a  background  distribution  with  B  0  can  be  generated 
by  adding  the  desired  \  alue  of  B  to  each  element  of  B,y. 

C.  BACKGROUND  MODEL  OUTPUTS 

The  background  model  was  developed  in  FORTRAN  code,  to  be  run  on  the  NTS 
IBM  3033  4381  mainframe  computer  network.  The  FORTRAN  code  for  this  model  is 
attached  as  Appendix  C. 

For  typical  IR  background  scenes  in  the  2-5fim  region,  as  determined  by  Itakura 
et  al.  [Ref  6],  a"'2s  170.0  mrad  and  a  varies  as  a  function  of  cloud  content.  Figure  2  on 
page  12  is  a  plot  of  the  amplitude  distribution  of  the  radiance  of  a  simulated  two- 
dimensional  background  generated  with  these  statistics,  and  covering  an  area  of  ap¬ 
proximately  2.1  mrad  by  418.9  mrad.  Likewise,  Figure  3  on  page  13  is  a  plot  of  the 
amplitude  distribution  of  the  radiance  of  a  different  simulated  two-dimensional  back¬ 
ground  covering  the  same  area,  but  generated  with  a  different  random  number  generator 
seed  and  a  dilferent  amplitude  variance. 

Both  of  these  backgrounds  have  been  chosen  to  be  acceptable  reperesentations  of 
infrared  backgrounds  in  the  7>-5nm  region,  because  their  amplitude  distributions  are 
similar  to  the  experimentally  determined  amplitude  distributions  shown  in  Figure  1  on 
page  7.  The  similarity  was  determined  on  the  basis  that  the  amplitude  distributions  of 
the  radiance  of  the  simulated  backgrounds  generally  conformed  to  the  shape  of  the  ex¬ 
perimentally  determined  amplitude  distribution  of  the  radiance  of  an  IR  background  in 
the  3-5nm  region.  This  can  be  seen  by  comparing  the  solid  and  dotted  plots  in 
Figure  2  on  page  12  and  Figure  3  on  page  13.  These  two  backgrounds  w^ill  be  used  in 
the  following  chapters  to  test  the  elTectiveness  of  the  LMS  filters. 

It  is  important  to  note  that  not  all  random  number  generator  seeds  produce  ac¬ 
ceptable  backgrounds.  In  fact,  approximately  50  different  seeds  were  tested  before  the 
two  backgrounds  discussed  above  were  generated.  Since  real  backgrounds  in  the 
3-5/im  ,  as  measured  by  Itakura  eial.  [Ref.  6],  typically  display  the  statistics  illustrated 
in  Figure  1  on  page  7,  simulated  backgrounds  not  displaying  a  basic  similarity  to  these 
statistics  were  considered  unacceptable  for  the  purpose  of  this  thesis. 
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OCCURRENCES  (NORMBLIZED) 


Figure  2.  Amplitude  Distribution;  Radiance  of  Simulated  Background  I:  Solid 
line,  amplitude  distribution  of  a  simulated  background  generated  with 
j(-‘=  170.0  mrad.  «r  =0.1,  and  random  number  seed  =356782367.  Dot¬ 
ted  line,  experimentally  determined  amplitude  distribution  of  an  IR 
background  in  the  3-5,um  region  (Ref.  6]. 
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Q  LD 


RMPLITUDE 


Figure  3.  Amplitude  Distribution;  Radiance  of  Simulated  Background  II:  Solid 
line,  amplitude  distribution  of  a  simulated  background  generated  with 
a-‘=  170.0  mrad,  a  =0.15,  and  random  number  seed® 41-1813567. 
Dotted  line,  experimentally  determined  amplitude  distribution  of  an  IR 
background  in  the  3-5A£m  region  (Ref  6]. 
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IV.  THE  IRSTD  MODEL 


A.  INTRODUCTION 

A  number  of  analytic  and  computer  simulated  performance  models  have  been  de¬ 
veloped  for  scanning  electro-optical  systems.  These  models  have  been  used  for  predict¬ 
ing  sensor  performance,  analyzing  system  transfer  functions,  and  developing  spatial 
filters. 

The  purpose  of  this  chapter  is  to  develop  a  computer  simulated  model  which  can 
be  used  to  study  the  NFS  Monterey  IRSTD  system.  The  model  closely  follows  the  an¬ 
alytic  model  presented  by  Scribner  and  Peters  (Ref.  9]. 

In  order  to  apply  linear-system  theory  to  this  model,  several  constraints  and  as¬ 
sumptions  are  necessaiy.  Specifically. 

•  the  system  must  be  linear  through  all  optical  processes, 

•  the  system  performance  must  be  representable  by  an  incoherent  response  function. 

•  the  optical  noise  sources  must  be  photon-additive, 

•  the  imaging  process  must  be  spatialy  invariant, 

•  the  system  bandwidth  must  be  narrow  enough  to  allow  for  reducing  spectrally  de¬ 
pendent  parameters  to  constants,  and 

•  the  background  clutter  must  be  stationary  and  Gaussian. 

The  simulation  models  one  detector  element,  scanning  horizontally  across  a  rectan¬ 
gular  object  plane.  .Although  this  model  represents  only  a  small  portion  of  the  NFS 
IRSTD  system,  it  can  still  be  used  to  test  and  develop  an  L.MS  filter  for  enhancing  the 
system  signal-to-noise  ratio. 

B.  FUNDAMENTAL  COMPONENTS  OF  THE  MODEL 

1.  The  Object  Plane 

The  object  plane  is  considered  to  be  a  two-dimensional  radiance  distribution, 
located  an  arbitrarj-  distance  from  the  IRSTD  aperture.  The  radiance  distribution  con¬ 
tains  information  from  both  targets  and  background  clutter,  and  can  be  represented  by 
the  function 


0(x',y')  =  T(x',y')-hB(x’,y') 


(4-1) 
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where  T(.\'.  y')  and  B(x'.  y')  describe  the  target  and  background  radiance  distributions 
respectively.  To  simplify  calculations  with  the  object  plane  distribution,  it  is  convenient 
to  define  an  ideal  object  plane  distribution  function  as 

0,(x,y)  =  T,(x,y)  +  B,(x,y)  (4  -  2) 

which  describes  the  object  plane  distribution  in  the  scale  of  the  image  plane. 

The  object  plane  power  spectral  density  (PSD)  is  simply  the  2-D  Fourier  trans¬ 
form  of  the  ideal  object  plane  distribution,  and  can  be  written  as 


0,(k,.  k,)  =  F2p{0,(x.y)} 


(4-3) 


2.  The  Image  Plane 

The  object  plane  radiance  distribution  is  mapped  onto  the  image  plane  by  the 
IRSTD  system.  This  process  can  be  represented  as  a  2-D  convolution  of  the  ideal  ob¬ 
ject  plane  radiance  distribution  function  0,(x.y)  .  with  the  system  point  spread  function 
P(\,y)  .  1  he  image  plane  radiance  distribution  function  can  therefore  be  written  as 

I(x,y)  =  0,(x.y)**Ps(x.y)  (4  -  4) 

wiiere  represents  the  2-D  convolution.  Additionally,  by  using  the  Convolution  Tlie- 
orem,  the  image  plane  PSD  can  be  written  as 


I(k,.  k  j  =  /-^^.illx.y)}  =  6,(k,.  kj  ) .  P^lk^,  kj)  (4  -  5) 

ITom  anah  sis  of  the  results  of  previous  studies  of  the  optical  performance 
characteri'>tics  of  the  NPS  IRSTD  ssstem  (Ref  lo],  it  can  be  assumed  that  the  system 
point  spread  function  can  be  represented  by  the  2-D  Gaussian  function 

P^(x.y)  =  — ?-y  y')  (4  -  6) 

2;tOQ 

where  is  the  radius  of  the  spot  size  of  a  point  target  in  the  focal  plane.  The  system 
response  function  is  the  Fourier  transform  of  the  point  spread  function,  and  can  be 
shown  to  be 


(4-7) 
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3.  The  1-D  Detector  Output 

The  detector  can  initially  be  considered  to  have  a  uniform  response  over  a  rec¬ 
tangular  surhicc.  with  the  dimensions  a  and  b  in  the  x  and  y  directions  respectively.  The 
detector  scans  the  image  plane  over  the  continuum  of  values  x  on  a  scan  line  centered 
at  >  =  b  2.  The  detector  output  can  then  be  expressed  as 


I  ^  2  )  fl*'*’  2  ...  , 

> )  =  I(x.y)dxdy 

T-t)  •’(5-f ) 


(4-8) 


Taking  the  limits  of  integration  over  x  into  the  integrand,  we  can  write 


Ic(x.  y) 


r(^'+f ) 

f^ 

l(x.y)  •  D(x  -  x.y)dx 

T-f) 

dy 


(4  -  91 


where  D(\  —  x.y)  is  the  detector  response  function. 

Since  the  inner  integral  in  equation  (4-9)  is  a  1-D  convolution  in  the  x-dircction. 
we  can  write,  using  to  represent  the  I-D  convolution. 


id(^-y)  = 


v) 

‘  'l(x,  y)«D(x.  y)dy 


(4  -  in) 


where  the  continuunt  of  values  \  is  identical  to  the  continuum  of  values  x.  and  tlie  de¬ 
tector  response  function  can  be  described  by 


Dix.y)  = 


Mi)  •«=«(!) 


(4-111 


If  the  detector  scans  the  image  plane  at  a  uniform  speed  the  output  can  be 
expressed  in  terms  of  time,  where  x  =  /?t.  Additionally,  if  the  detector  responsivity  is 
known,  the  detector  output  voltage  can  be  written  as 

V(t)  =  5?.M/lt,y)  (4-12) 
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The  detector  output  voltage  can  also  be  expressed  in  the  electrical  frequency 
domain  as 


V{0  =  F,^{^.Ia(/?t,y)} 


(4-13) 


4.  Noise  and  AC  Coupling 

One  type  of  noise  of  particular  interest  in  the  IRSTD  system  is  white  noise. 
Considering  this  noise  to  be  additive,  the  detector  output  can  be  expressed  as 

V'(t)  =  V(t)  +  NJt)  (4-14) 

where  N^(t)  describes  the  white  noise,  which  can  be  represented  by  a  uniformly  distrib¬ 
uted  random  variable. 

Since  the  outputs  are  AC  coupled,  it  is  of  interest  to  describe  this  coupling  in 
the  model.  The  AC  coupled  outputs  can  be  expressed  as  [Ref  10) 

Vc(t)  =  V'(t)  V'(t)"e-‘ ''''  (4-15) 

where  Qo  is  the  capacitor  charge  at  t  =  0,  R  is  the  resistance,  C  is  the  capacitance  of  the 
AC  coupling  circuit,  and  *  again  represents  a  1-D  convolution.  As  a  first  approxi¬ 
mation,  Qo  can  be  considered  to  be  equal  to  zero;  then 

Vc(t)  =  V'(t)  -  V'(t)^e-‘ ''''  (4-16) 

approximately  describes  the  output  of  the  detector. 

.At  this  point,  equation  (4-16)  describes  the  output  of  a  generic  scanning  IR 
sensor  system.  By  including  parameters  representative  of  the  NFS  IRSTD  system,  this 
model  can  be  used  to  simulate  IRSTD  detector  output  data.  These  specific  parameters 
will  be  introduced  in  the  following  section. 

C.  IRSTD  MODEL  OUTPUTS 

This  model  was  also  developed  in  FORTRAN  code,  and  run  on  the  NFS  IBM 
3033  4381  computer.  The  FORTRAN  code  is  attached  as  Appendix  D.  The  program 
was  developed  to  produce  IRSTD  output  signals  ranging  in  amplitude  from  -0.35  to 
0.35.  This  was  done  because  digital  data  from  a  12-bit  analog-to-digital  (A/D)  converter 
can  be  easily  used  to  represent  output  signals  ranging  in  amplitude  from  -0.4095  to 
0.4095.  Specific  parameters  used  in  the  simulation  include:  detector  width  a  =  0.3  mrad, 
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detector  height  b=  2.0  mrad,  spot  size  radius  =  0.7  mrad,  and  AC  coupling  constant 
RC=  1.3.  These  specific  parameters  were  used  to  produce  output  data  representative 
of  real  IRSTD  output  data.  It  must  be  emphasized  that  the  model  was  not  designed  to, 
nor  was  it  intended  to,  produce  an  "exact"  replica  of  real  IRSTD  data. 

Insening  backgrounds  1  and  11  into  this  IRSTD  model,  results  in  the  detector  out¬ 
puts  illustrated  in  Figure  4  on  page  19  and  Figure  5  on  page  19.  These  two  detector 
outputs  will  be  used  throughout  this  thesis,  and  will  be  referred  to  as  BKDl  and  BlCDll 
respectively. 

For  purposes  of  comparison.  Figure  6  on  page  20  is  a  real  IRSTD  background  ob¬ 
tained  from  the  NFS  IRSTD  system  during  sj-stem  calibration,  (Ref.  2].  This  back¬ 
ground  contains  information  from  a  distant,  uniform,  partly-cloudy  sky,  with  ambient 
atmospheric  temperature  of  52.5  degrees  F  and  visibility  of  approximately  20  miles. 
Additionally,  the  data  vvas  digitized  at  100  ^s  per  sample  on  a  Nicolet  digital 
oscilloscope,  with  the  system  scanning  at  n  radians  per  second.  This  background  will 
hereupon  be  referred  to  as  BKDIII. 

In  addition  to  generating  simulated  IRSTD  outputs  of  background  distributions, 
this  model  was  also  used  to  generate  simulated  IRSTD  outputs  of  various  target  signals. 
These  simulated  target  signals  were  generated  by  inserting  target  signals  of  uniform 
amplitude  and  finite  dimensions  into  an  ideal  object  plane  distribution.  Equation  (4-2), 
having  a  background  distribution  identically  equal  to  zero. 

To  generate  IRSTD  outputs  simulating  point  and  near  point  targets,  target  signals 
of  0.10472  mrad  by  0.10472  mrad,  0.31416  mrad  by  0.3I4I6  mrad,  0.52360  mrad  by 
0.52360  mrad,  1.0472  mrad  by  1.0472  mrad,  and  1.5708  mrad  by  1.5708  mrad  were  used. 
The  simulated  1 RSTD  detector  outputs  of  these  five  targets  are  illustrated  in  Figure  7 
on  page  20,  Figure  8  on  page  21,  Figure  9  on  page  21,  Figure  10  on  page  22,  and  Fig¬ 
ure  1 1  on  page  22, 

To  show  how  well  the  model  represents  the  IRSTD  system,  these  model  detector 
outputs  of  simulated  targets  can  be  compared  to  the  detector  output  of  a  whip  antenna 
taken  by  the  NFS  IRSTD  system,  Figure  12  on  page  23.  These  five  simulated  targets 
and  one  real  target  will  be  used  throughout  this  thesis,  and  will  be  referred  to  as  TGTI, 
TGTII,  TGTIII,  TGTIV,  TGTV,  and  TGTVI  respectively. 
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Figure  4.  Detector  Output;  Simulated  Background  #  1  (BKDl) 
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Figure  6.  Detector  Output;  Real  Background  (BKDIII) 
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Figure  7.  Detector  Output;  Simulated  Target  ^1  (TGTI):  Target  dimensions  arc 
0.10472  mrad  by  0.10472  mrad. 


20 


Figure  9.  Detector  Output;  Simulated  Target  #3  (TGTIII):  Target  dimensions 
are  0.52360  mrad  by  0.52360  mrad. 
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Figure  10.  Detector  Output;  Simulated  Target  #-4  (TGTIV):  Target  dimensions 
are  1.0472  mrad  by  1.0472  mrad. 


Figure  11.  Detector  Output;  Simulated  Target  #5  (TGTV'):  Target  dimensions 
are  1.5708  mrad  by  1.5708  mrad. 
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_ nZIMUTH  (MRnP) _ 

Figure  12.  Detector  Output;  Real  Target  (Whip  Antenna)  (TGTVT) 
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V.  THE  LMS  FILTER 


A.  THEORY 

As  previously  discussed,  the  spatial  extent  of  background  clutter  sources  is  typically 
much  greater  than  the  spatial  extent  of  targets  (Ref  3).  This  can  be  shown  by  comparing 
the  frequency  content  of  the  spatial  frequency  distributions  of  the  backgrounds  and 
targets.  Figure  13  through  Figure  17  on  page  26  are  plots  of  the  magnitudes  of  the 
spatial  frequency  distributions  of  BKDI,  BKDII,  TGTI,  TGTV,  and  TGTVT  respec¬ 
tively.  These  figures  clearly  illustrate  that  both  the  simulated  and  real  background  dis¬ 
tributions  are  of  predominantly  lower  frequency  content,  and  hence  larger  spatial  extent, 
than  the  target  signals. 


o 

a 


Figure  13.  Magnitudes  of  the  Spatial  Frequency  Distribution  of  BKDI 
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Figure  16.  Magnitudes  of  the  Spatial  Frequency  Distribution  of  TGTV 


Figure  17.  Magnitudes  of  the  Spatial  Frequency  Distribution  of  TGTVI 
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A  filter  which  has  been  shown  to  be  very  effective  under  these  conditions,  can  be 
developed  using  the  principle  of  least-mcan-square  (LMS)  optimization  [Ref.  3].  The 
response  of  such  a  filter  can  be  determined  by  minimizing  the  difference  between  the 
received  signal  and  a  function  which  approximates  the  presence  of  both  target  and 
background.  The  most  crucial  step  in  this  technique  is  selecting  the  proper  function  for 
estimating  the  target  signal  and  background. 

The  model  signal  which  best  estimates  the  presence  of  both  target  and  background 
clutter  can  be  expressed  as  a  weighted  sum  of  anticipated  target  and  an  expansion  series 
of  low  frequency  noise  expanded  about  the  center  of  the  target  signal 

x(t)  =  a  s(t)  +  bo  +  b,T  +  biT^  +  •  •  •  (5-1) 

where  s(t)  is  the  model  target  signal,  t  is  time  measured  from  the  center  of  s(t),  and  a, 
bij.  b,.  b,.  .  .  .  are  weighting  coefficients  specifying  the  amount  each  component  of  the 
series  contributes  to  the  signal. 

Defining  the  function  v(t)  as  an  actual  signal  output  from  the  IRSTD  system,  the 
presence  of  a  target  at  a  given  test  time  to  can  be  tested  for  by  minimizing  the  squared 
sum 


F(to)  =  I '  dt[v(to  -  t)  -  x(t)]‘  (5  -  2) 

where  t^  is  the  half-amplitude  width  of  the  model  target  signal  s(t). 

Since  the  coefficient  a  .  in  Equation  (5-1),  describes  the  amount  of  target  present  in 
v(t)  at  t  =  t(,.  it  can  be  considered  to  be  the  time  dependent  output  of  an  optimum  target 
detection  LMS  filter.  To  simplify  the  minimization  of  equation  (5-2),  it  is  convenient  to 
truncate  equation  (5-1 ).  If  we  consider  removing  only  DC  offset  and  linear  noise  from 
the  recei%ed  signal,  equation  (5-2)  can  be  rewritten  as 

F(to)  =  [  '  dt[v(to  -  t)  -  a  s(t)  -  bo  -  b,T]^  (5  -  3) 

•’-td 

Differentiating  equation  (5-3)  with  respect  to  the  coefficients  a,  bo ,  and  b„  and  set¬ 
ting  each  derivative  equal  to  zero,  leads  to  the  equations 

r}' 

— : —  =  <  sv  >  —  a  <  s'  >  —  bn  <  s  >  —  bt  <  ST  >  =  0  (5  —  4) 

f  a 
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cF  2 

-t; —  =  <v>-a<s>-bo-b,  <T  >  =  0 
cbo 


rF  2 

=  <vt>-a<ST>-bo<T>-b,  <T  >  =  0 

C  Di 


(5-5) 

(5-6) 


where  <  >  denotes  the  time  average  over  the  interval  —  t^  t  ^  t^.  Equations  (5-4), 
(5-5),  and  (5-6)  can  be  written  in  matrix  form  as 


<s>  <s>  <ST>  a 

<  s  >  1  <  T  >  bo 

<  ST  >  <  T  >  <  T^  >  b, 


<  sv  > 

<  V  > 

<  VT  > 


(5-7) 


or,  since  <  st  >  =  0  and  <  t  >  =  0  for  symmetric  target  signals 

r_i  _  A  1  r _ 


a  <  s  >  <  s  > 
bo  =  <  s  >  1 

b,  0  0 


0  <t"> 


<  SV  > 

<  V  > 


(5-8) 


Then,  since  the  inverse  of  the  3x3  matrix  in  Equation  (5-S)  can  be  described  by  the 
proportionality 


<  s‘  > 

<  s  > 

0 

<  > 

—  <  S  >  <  > 

0 

<  s  > 

1 

0 

-'oc  -<S><t'> 

<  S^  ><  T^  > 

0 

0 

0 

A 

V 

1 _ 

0 

0 

2  2 
<  s  >  —  <  s  > 

Equation  (5-8)  can  be  rewritten  as 


<  T  > 


b„  OC  -  <  S  ><  T^  >  <  S^  ><  > 


0  <  sv  > 

0  <  V  > 

<  s^  >  -  <  S  <  VT  > 


(5-9) 


Solving  for  a  results  in  the  expression 


a  OC  <T^><SV>  —  <T^><S><V> 


(5  -  10) 
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Rewriting  the  time  a\  erages  in  integral  form,  and  recognizing  that  dividing  the  right  side 
of  Equation  (5-10)  by  <t‘>  will  not  destroy  the  proportionality,  leads  to 


td  td  f‘d  fq 

—  J  dt  s(t)  v(to  - 1)  -  —  J  dt  s(t)  J  dt  v(to  - 1) 


This  expression  can  be  simplified  to 


~J  '  dt[s(t)  -  <  s  >  ]v(to  -  t) 


a  o:  -J-J  "  dt  h(t)  v(to  -  t) 


(5-11) 


where  h(t)  is  the  LMS  filter  function 


h(t)  =  s(t)  -  <  s  > 


(5-12) 


B.  FILTER  DEVELOPMENT 

In  this  section,  an  LMS  filter  will  be  developed  for  extracting  target  signals  from 
IRSTD  outputs.  This  is  accomplished  by  properly  modeling  IRSTD  target  signals,  and 
generating  the  filter  function  from  equation  (5-12). 

.As  previously  mentioned,  the  most  crucial  step  in  this  filtering  technique  is  choosing 
the  proper  function  to  model  target  signals.  In  order  to  use  equation  (5-12)  to  develop 
this  filter,  the  target  signal  must  be  modeled  as  a  symmetric  function.  Referring  back 
to  the  target  signals  illustrated  in  Chapter  IV,  it  is  clear  that  these  signals  are  not  sym¬ 
metric.  However,  upon  closer  inspection,  it  can  be  seen  that  each  target  signal  is  ap¬ 
proximately  symmetric  over  a  small  region  about  its  maximum  amplitude. 

Four  functions  which  have  the  potential  of  modeling  this  symmetric  region  of  the 

target  signal  are:  a  raised  cosine  function,  a  parabolic  function,  a  Gaussian  function, 

^  sin(7rxi  .  .  .  .  ,  , 

and  a  Sine  function  (Sinc(x)  =  — ; — — ).  Upon  initial  investigation,  the  parabolic, 

Gaussian,  and  Sine  functions  were  ruled  out,  because  the  parabolic  function  w’ould  not 

fit  the  top  region  of  the  target  signals,  and  the  Gaussian  and  Sine  functions  would  not 
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fit  the  sides  of  the  target  signals.  The  raised  cosine  function,  however,  modeled  the  tar¬ 
get  signals  remarkably  well,  over  a  range  near  the  peak  of  the  signals,  and  was  selected 
as  an  appropriate  function  for  modeling  them. 

Figure  IS  on  page  31,  through  Figure  25  on  page  34  display  how  raised  cosine 
functions  of  half-amplitude  widths  ranging  from  0.31416  to  1.04728  mrad  can  be  used 
to  model  TGTl.  It  can  be  seen  from  these  figures  that  the  raised  cosine  functions  of 
half-amplitude  widths  ranging  from  0.31416  to  0.94248  mrad  model  the  leading  edge,  or 
rise,  of  the  target  signal  remarkably  well.  However,  raised  cosine  functions  with  half¬ 
amplitude  widths  greater  than  0.62832  mrad  begin  to  deviate  substantially  from  the 
trailing  edge,  or  fall,  of  the  target  signal. 

Similarly,  upon  investigating  how  well  raised  cosine  functions  model  TGTll  through 
TGTVl.  it  was  found  that  raised  cosine  functions  of  half-amplitude  widths  ranging  from 
0.31416  to  0.94248  mrad  also  modeled  the  leading  edges  of  these  target  signals  remark¬ 
ably  well.  However,  raised  cosine  functions  with  half-amplitude  widths  greater  than 
0.62S32  nuad  begin  to  deviate  substantially  from  the  trailing  edge  of  TGTll  and 
TGI  111.  functions  with  half-amplitude  widths  greater  than  0.73304  mrad  begin  to  devi¬ 
ate  substantially  from  the  trailing  edge  of  TGTIV.  and  functions  with  half-amplitude 
widths  greater  than  0.83776  mrad  begin  to  deviate  substantially  from  the  trailing  edge 
of  TGTV  and  TGTVl. 

From  this  analysis,  it  can  be  concluded  that  the  symmetric  portion  of  the  target 
signals  can  be  modeled  as  a  raised  cosine  function 

=  0  elsewhere 

where  t^  is  the  half-amplitude  width  of  the  signal,  also  refered  to  as  the  signal  dwell  time. 

It  should  be  noted  that  t^  can  be  expressed  in  terms  of  either  time  or  azimuth,  l.'n- 
less  otherwise  specified,  it  will  always  be  expressed  in  terms  of  azimuth,  with  units  of 
mrad.  throughout  the  remainder  of  this  thesis. 


30 


rmplitude:  ^  rmplitude 


8  18.  TGTI  and  Raised  Cosine  with  0.3I4I6  nirad  Half-Amplitude  Width 


LO 


Figure  19.  TGTI  and  Raised  Cosine  with  0.4 1888  mrad  Half- Amplitude  Width 
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Figure  20.  TGTI  and  Raised  Cosine  with  0.52360  mrad  Half- Amplitude  Width 


Figure  21.  TGTI  and  Raised  Cosine  with  0.62832  mrad  Half- Amplitude  Width 


Figure  22.  TGTI  and  Raised  Cosine  with  0.73304  mrad  Half-Amplitude  Width 


Figure  23.  TGTI  and  Raised  Cosine  with  0.83776  mrad  Half-Amplitude  Width 
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Figure  24.  TGTI  and  Raised  Cosine  with  0.94248  mrad  Half- Amplitude  Width 
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Figure  25.  TGTI  and  Raised  Cosine  with  1.0472  mrad  Half-Amplitude  Width 


Substituting  equation  (5-13)  into  equation  (5-12)  gives 


h(t)  =  cos7r(  -^ )  —  tj  <  t  ^  td 

td  ^  1  (5-14) 

=  0  elsewhere 

where  <  1  +  cosn{  -p )  >  =  1- 

Since  the  output  signals  of  the  IRSTD  system  are  digitized,  it  is  necessary'  to  express 
the  filter  in  digital  form.  This  can  be  accomplished  by  expressing  s(t)  in  the  form  of  a 
vector,  whose  components  are 


in  =  y  1  +  cos7m(  ) 


(5-15) 


where 


n  =  0,  i  l*i2, •••  +  ” 


(5-16) 


and  t,  is  the  sampling  time  of  the  IRSTD.  It  should  be  noted  that,  if  -f-  is  non-integer, 
the  maximum  value  of  n  will  assume  the  value  of  the  integer  which  is  closest  to  and  less 
than  For  example,  if  -^  =  3.2;  n  =  0.  ±1,  ±2,  ±3. 

To  illustrate  this  vector  representation  for  —  =  3,  it  is  clear  that 

Sn  =  ^[l -Hcosrr-^]  n  =  (-3, -2, -1,  0,  1,  2,3)  (5-17) 


or  in  matrix  form 


S  =  — 
^  2 


0.000 

0.000 

0.500 

0.250 

1.500 

0.750 

2.000 

1.000 

1.500 

0.750 

0.500 

0.250 

0.000 

0.000 

(5-18) 


The  vector  representation  of  the  filter  can  then  be  written,  from  equation  (5-12),  as 

1 


H  = 


|Smax-<S> 


[S  -  <  S  >  ] 


(5-19) 
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where  the  factor  -r^ - - — r  is  a  normalization  constant.  To  continue  the  example 

t,  |S^„-<S>| 

for  7“  =  3,  the  LMS  filter  coefiicients  become 


-0.429 

-0.750 

-0.179 

-0.312 

0.321 

0.562 

0.571 

= 

1.000 

0.321 

0.562 

-0.179 

-0.312 

-0.429 

-0.750 

(5  -  20) 


This  example  LMS  filter  H  is  symmetric  with  zero  mean,  making  it  appropriate  for 
extracting  symmetric  target  signals  from  received  signals  digitized  with  a  sampling  time 
of  t,  =  3td,  and  containing  linear  noise  of  spatial  extent  much  greater  than  t^. 

When  the  limits  of  integration  are  changed  to  investigate  the  presence  of  target 
signals  over  all  time,  the  coelTicient  a  becomes  time  dependent,  and  equation  (5-12)  can 
be  rewritten  as 


a(to) 


td  f"' 

"tJ 


dt  h(t)  v(to  - 1) 


(5-21) 


which  is  the  convolution  of  v(t)  with  h(t).  In  matrix  form,  equation  (5-21)  can  be  re¬ 
written  as 


A^H‘V  (5  -  22) 

Therefore,  the  output  array.  A,  can  be  determined  by  digitally  convolving  the 
IRSTd  output  signal.  V,  with  the  digital  LMS  filter,  H.  It  should  be  emphasized  that 
the  output  array.  A,  is  the  vector,  or  digital,  representation  of  the  time  dependent  ex¬ 
pansion  series  coefficient  a,  which  describes  the  amount  of  target  present  in  V. 

For  large  IRSTD  output  arrays,  V,  digital  convolution  can  be  costly  with  respect  to 
computer  processing  time.  It  is  therefore  convenient  to  take  advantage  of  the  Convo¬ 
lution  Theorem,  and  rewrite  equation  (5-22)  as 

A=F7i{H.V}  (5-23) 
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where  H  and  V  are  the  Fourier  transforms  of  H  and  V  respectively,  Fji  represents  the 
ID  inverse  Fourier  transform,  and  H  is  commonly  referred  to  as  the  filter  transfer  func¬ 
tion. 

In  order  to  accomplish  the  spatial  frequency  domain  multiplication  of  equation 
(5-23).  it  is  first  necessar>-  to  place  the  coefficients  of  the  LMS  filter  in  an  array  of  the 
same  dimensions  as  V.  This  is  done  by  centering  the  coefficients  in  an  empty  array,  and 
letting  all  of  the  remaining  elements  of  the  array  be  zero. 

A  set  of  trial  filters,  for  filtering  both  simulated  and  real  IRSTD  data,  will  be  gen¬ 
erated  in  the  following  section.  The  resultant  filtered  outputs  will  be  presented  and  dis¬ 
cussed  in  the  following  chapters. 

C.  TRIAL  FILTERS 

As  previously  discussed,  a  raised  cosine  function  can  be  considered  an  appropriate 
model  for  specific  portions  of  TGTI  through  TGTVT.  To  see  how  well  this  function 
does  model  the  targets,  a  set  of  15  trial  LMS  filters  were  developed  from  raised  cosine 
functions  with  half-amplitude  widths  ranging  from  0.10472  to  1.5708  mrad. 

Using  the  methods  outlined  in  the  previous  section,  the  coefficients  of  these  trial 
filters  were  determined,  and  were  used  to  filter  both  simulated  and  real  IRSTD  output 
data  with  the  objective  of  seeing  which  filters  yield  the  best  signal-to-noise  (S  \)  ratio 
increase  for  each  target.  It  is  important  to  note  that  even  though  these  filters  were  de¬ 
signed  for  extracting  targets  with  specific  signal  half-amplitude  widths,  they  are  much 
more  general  in  nature.  Since  the  targets  under  investigation  were  sampled  at  0.10472 
mrad  per  sample,  these  filters  are  actually  capable  of  extracting  any  targets  with  signals 
having  1  through  15  samples  per  half  amplitude  width.  In  light  of  this,  these  filters  will 
hereupon  be  referred  to  as  LMSl  through  LMS15  respectively. 

The  coefficients  of  LMSl  through  LMS8  are  listed  in  Table  1  on  page  39,  and  the 
coefficients  of  LMS9  through  LMSl 5  are  listed  in  Table  2  on  page  40.  The  coefficients 
of  a  particular  filter  can  be  obtained  from  these  tables  by  reading  down  the  columns 
under  the  appropriate  heading  for  the  desired  filter.  The  coefficients  were  rounded-off 
to  the  third  significant  digit  to  the  right  of  the  decimal  for  ease  in  illustration.  However, 
the  coefficients  used  in  the  filtering  process  were  generated  by  the  NFS  IBM  3033  4381 
computer  and  represented  by  32-bit  words,  thereby  increasing  their  precision  substan¬ 
tially. 
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Additionally,  it  should  be  noted  that  by  using  these  LMS  filters  on  TGTI  through 
TGTVI.  an  inherent  thresholding  takes  place  during  the  filtering  process.  This  thresh¬ 
olding  can  be  seen  by  referring  back  to  Figure  18  on  page  31  through  Figure  25  on  page 
34,  where  the  target  signals  can  be  considered  to  be  thresholded  at  the  base  of  the  raised 
cosine  functions  during  the  filtering  process.  Therefore,  as  the  filter  width  increases,  a 
greater  portion  of  target  signal  becomes  available  for  filtering.  This  combination  of  fil¬ 
tering  and  thresholding,  during  the  filtering  process,  can  generate  serious  problems  in 
that  the  probability  of  thresholding  out  target  signals  increases  as  the  width  of  the  raised 
cosine  function  decreases.  It  will  therefore  be  necessary'  to  weigh  this  factor  when  de¬ 
termining  which  filter  is  optimal,  by  considering  both  the  noise  reduction  capability  and 
half-amplitude  width  of  the  filters. 

Before  proceding  to  the  next  chapter,  a  discussion  of  the  real  background,  BKDIII, 
is  necessary  .  As  previously  mentioned,  both  the  IRSTD  model  and  the  IRSTD  DCL 
sample  data  at  0.10472  mrad  per  sample.  However,  the  real  background  was  sampled 
by  the  Xicolet  digital  oscilloscope  at  0.31416  mrad  per  sample,  or  one  third  the  system 
sampling  rate. 

In  order  to  embed  the  targets  into  BKDIII,  it  was  necessary  to  sample  them  at  the 
same  rate  as  the  background.  This  was  accomplished  by  altering  the  sampling  rate  of 
the  IRSTD  model,  and  generating  new  targets  sampled  at  0.31416  mrad  per  sample. 
.•\nalysis  of  these  modified  target  signals  showed  that  the  same  raised  cosine  functions 
used  to  model  the  simulated  targets  can  be  used  to  model  the  modified  targets  as  well. 
However,  since  the  increased  sampling  time  reduces  -p  by  a  factor  of  three,  the  LMS 
filters  will  not  have  the  same  effect  on  BKDIII  as  they  do  on  BKDI  and  BKDII.  Spe¬ 
cifically,  if  LMS9  is  determined  to  be  the  optimal  filter  for  TGTI  in  BKDI  and  BKDII, 
L.MS3  can  be  e.xpected  to  be  optimal  for  TGTI  in  BKDIII. 
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Table  1.  LMS  FILTER  COEFFICIENTS _ 

Azimuth _ LMS  Filter 


(mrad) 

LMSl 

LMS2 

LMS3 

LMS4 

LMS5 

LMS6 

LMS7 

LMS8 

mnm 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

-0.889 

Is 

U.OtX) 

0.000 

o.ooo 

0.000 

0.000 

0.000 

-0.875 

-0.817 

-6  [, 

o.ooo 

0.000 

0.000 

0.000 

0.000 

-0.857 

-0.782 

-0.612 

Is 

O.U<  HI 

0.000 

0.000 

0.(X)0 

-0.833 

-0.733 

-0.522 

-0.306 

-4  r 

O.IHIO 

o.ooo 

0.0(XI 

-o.soo 

-0.658 

-0.393 

-0. 146 

0.056 

Is 

0.(HHI 

0.000 

-0.750 

-0.536 

-0.200 

0.071 

0.271 

0.417 

■-  ^ 

0.000 

-0.667 

-0.312 

0.100 

0.367 

0.536 

0.647 

0.723 

-1  Is 

-0.500 

0.167 

0'36 

0.825 

0.876 

0.907 

0.928 

0  !. 

1.000 

1.000 

1.0o(.) 

l.(XX> 

1.000 

1 .000 

1 .000 

1  .(XX_) 

1 1. 

-0.500 

0.167 

0.736 

0.825 

0.876 

0.907 

0.928 

:  /. 

o.ni  10 

-0.667 

-0.312 

0.100 

0.367 

0.536 

0.647 

0.723 

O.ooo 

O.I.H.X) 

-0.750 

-0.536 

-0.200 

0.071 

0.271 

0.417 

4  I. 

O.IM.lO 

O.OIHI 

0.000 

-0.8(X) 

-0.65S 

-0.146 

0.056 

>  ts 

0.(H  HI 

0.o0(i 

0.(H  10 

0.000 

-0.733 

-0.522 

-0.3(»6 

6  t. 

O.o(  Id 

0.000 

0.0(X> 

0.000 

0.000 

-0.857 

BIO 

-0.612 

7  /, 

O.ooo 

0.0(10 

0.(100 

0.(.Xi0 

O.O(X) 

O.O(X) 

-0.875 

-(J.817 

S  t. 

O.Ol  10 

0.(  XH.I 

0.000 

O.(XX) 

O.CMK) 

0.000 

0.()(.)0 

-O.SS9 

Table  2. 

LMS  FILTER  COEFFICIENTS 

Azimuth 

LMS  Filter 

(mrad) 

LMS9 

LMSIO 

LMS  11 

LMS12 

LMS  13 

LMS  14 

LMS15 

-15  r, 

*).(,)(«.» 

0.000 

0.000 

0.000 

0.000 

0.000 

-0.937 

-14  r 

0.000 

0.000 

0.000 

0.000 

0.000 

-0.933 

-0.916 

-13  i, 

O.oOO 

0.000 

0.000 

0.000 

-0.929 

-0.909 

-0.854 

■12  r, 

0.000 

0.000 

0.000 

-0.923 

-0.901 

-0.838 

-0.752 

-11  t, 

0.00(.) 

0:000 

-0.917 

-0.890 

-0.818 

-0.722 

-0.617 

■10 

0.000 

-0.909 

-0.878 

-0.794 

-0.686 

-0.569 

-0.453 

-9/, 

-0.900 

-0.S62 

-0.765 

-0.641 

-0.512 

-0.386 

-0.268 

-S  t. 

-0.843 

-0.727 

-0.586 

-0.442- 

-0.306 

-0.182 

-0.070 

■1  i. 

-U.67S 

-0.516 

-0.356 

-0.210 

-0.081 

0.033 

0.133 

■b  t. 

-0.425 

-(•.250 

-0.095 

0.038 

0.152 

0.248 

0.331 

-5  /. 

-0.115 

0.045 

0.178 

0.287 

0.378 

0.453 

0.516 

-4  [. 

(t.215 

0.34(» 

0.440 

0.519 

0.583 

0.636 

0.679 

-3  /. 

0.525 

0.607 

0.669 

0.718 

0.757 

0.789 

0.815 

-:  1. 

0.~78 

0.818 

0.848 

0.871 

0.890 

0.904 

0.916 

■  -1 1. 

0.943 

msm 

0.961 

0.967 

0.972 

0.976 

0.979 

0  / 

1.  (»(.«) 

l.tw)0 

1 .000 

1.000 

1.000 

1.0(,)0 

1.000 

1 1. 

BOBU 

0.953 

0.961 

0.967 

0.972 

0.976 

0.979 

2  t 

0.778 

0.818 

0.848 

0.871 

0.890 

0.904 

0.916 

3  t. 

0.525 

0.6(t7 

0.669 

0.718 

0.757 

0.789 

0.815 

b 

0.215 

0.340 

0.440 

0.519 

0.583 

0.636 

0.679 

5  [. 

-0.115 

0,(»45 

0.178 

0.287 

0.378 

0.453 

0.516 

b  1 

-0.425 

-0.250 

-0.095 

0.038 

0.152 

0.248 

0.331 

7  i. 

-0.67S 

-0.516 

-0.356 

-0.210 

-0.081 

0.033 

BBIHM 

b 

-0.727 

-0.586 

-0.442 

-0.3(.Xi 

-0.182 

-0.070 

9  i. 

.0.9(10 

-0.862 

-0.765 

-0.641 

-0.512 

-0.386 

-0.268 

1(»  i. 

0.000 

-0.909 

-0.878 

-0.794 

-0.686 

-0.569 

-0.453 

11  /, 

0.000 

0.000 

-0.917 

-0.890 

-0.818 

-0.722 

-0.617 

12  L 

O.OoO 

0.000 

0.000 

-0.923 

-0.901 

-0.838 

-0.752 

1 3  r 

0.(»i.>(,» 

0.000 

0.000 

0.000 

-0.929 

-0.909 

-0.854 

14  /, 

0.000 

0.0(>0 

0.000 

0.000 

0.000 

-0.933 

-0.916 

15  i. 

0. ()<»(,) 

O.oOO 

0.(XK) 

0.000 

0.000 

0.000 

-0.937 

VI.  FILTER  OUTPUTS 


The  15  trial  LMS  filters  (LMSl  through  LMS15),  introduced  and  developed  in  the 
previous  chapter,  were  used  to  filter  both  simulated  and  real  IRSTD  outputs.  The  fol¬ 
lowing  sections  summarize  the  effectiveness  of  these  filters,  and  illustrate  filtered  outputs 
of  selected  LMS  filters. 

A.  SIGNAL-TO-NOISE  RATIOS 

In  order  to  quantitatively  analyze  the  effectiveness  of  the  LMS  filters,  it  is  necessary 
to  determine  the  signal-to-noise  {S  N)  ratio  of  both  the  filtered  and  unfiltered  outputs. 
Since  the  convolution  of  the  filter  with  the  IRSTD  output  is  performed  in  the  spatial 
frequency  domain,  it  is  convenient  to  determine  the  S,  N  ratio  in  this  domain  also. 

The  S  N  ratio  of  the  filtered  detector  outputs  can  be  represented  by  the  expression 


(H..T.) 

(S  N)p  =  - -  (6-1) 

1  =  0 

where  T,  are  the  elements  of  the  target  spatial  frequency  distribution,  B,  are  the  elements 
of  the  background  spatial  frequency  distribution,  and  the  index  i  ranges  over  all  of  the 
N  elements  in  the  IRSTD  output  array. 

It  is  important  to  note  that  for  analysis  of  the  S  X  ratios  of  the  unfiltered  IRSTD 
data,  no  filter  is  being  used.  Therefore,  the  elements  of  the  filter  transfer  function,  H, 
are  all  equal  to  one.  and  Equation  (6-1)  can  be  rewritten  as 
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n 


(S.N)u - p -  (6-2) 


for  the  unfiltered  IRSTD  data. 

B.  ERROR  ANALYSIS 

In  order  to  obtain  some  degree  of  confidence  in  the  filtered  outputs,  it  is  essential 
to  study  the  types  of  errors,  or  uncertainty,  which  are  inherent  in  computer  simulations 
and  calculations.  Two  of  the  most  important  sources  of  uncertainty,  which  are  appli¬ 
cable  to  this  thesis,  are  quantization  and  roundoff. 

Quantization  error  occurs  during  A  D  conversion,  where  digital  quantities  can  only 
be  represented  to  a  fixed  and  finite  precision.  The  effect  of  this  limitation  to  finite 
quantities  is  that  the  representation  of  all  physical  and  simulated  signals  must  be  viewed 
as  corrupted  by  noise.  The  variance  of  this  noise  can  be  represented  by  the  expression 
[Ref  111 


(6  -  3) 


where  b  is  the  number  of  bits  being  used  to  represent  the  magnitude  of  the  digital  data. 

Roundoff  error  is  introduced  by  rounding  off  the  least  significant  figure  kept,  in  or¬ 
der  to  represent  digital  quatities  in  a  fixed  and  finite  form.  The  effect  of  this  error  must 
be  viewed  as  an  addition  of  a  noise  term,  with  a  variance  which  can  be  represented  by 
the  expression  [Ref  1 1] 


=  (6-4) 

where  d  is  the  number  of  base  ten  significant  digits  to  the  right  of  the  decimal  point  re¬ 
maining  after  roundoff. 

The  effects  of  quantization  and  roundoff  errors  on  the  filtered  outputs  can  be  de¬ 
termined  by  considering  the  noise  amplification  effects  of  digital  filters.  Specifically,  if 
the  variance  of  the  elements  of  the  unfiltered  signal  (aj)  is  known,  and  if  the  variance  in 
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the  individual  filter  coefficients  is  negligible  compared  to  trj,  the  variance  of  the  elements 
of  the  filtered  output  ((t;)  can  be  determined  from  [Ref  12] 


(6-5) 


where  H,  are  the  filter  coefficients,  and  the  sum  is  over  the  N  elements  of  the  IRSTD 
output  signal. 

For  the  purpose  of  this  thesis,  all  digital  elements  of  the  target  and  background 
distributions  were  either  quantized  or  simulated  to  be  quantized  by  a  13-bit  A/D  con¬ 
verter.  Of  these  13  bits,  one  bit  is  used  to  designate  sign  and  12  bits  are  used  to  specify 
magnitude.  Therefore,  the  variance  of  the  elements  of  the  background  distributions  and 
target  signals  can  be  considered  to  be  4.97  x  10'’.  Since  the  target  signals  are  imbedded 
into  the  background  distributions  to  construct  the  unfiltered  signal,  the  variance  of  the 
elements  of  the  unfiltered  signal  must  be  represented  as  the  sum  of  the  variances  of  the 
elements  of  the  background  distributions  and  target  signals,  which  is  equal  to 
9.94  X 

Using  Equation  (6-5),  the  variances  of  the  elements  of  the  filtered  outputs  for  LMSl 
through  IMS  15  were  determined,  and  are  listed  in  Table  3  on  page  44.  These  variances 
give  a  measure  of  the  noise  amplification  efTect  of  the  filtering  process,  which  becomes 
more  and  more  significant  as  the  width  of  the  filter  increases. 
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Table  3.  3'ARIANCES  OF  ELEMENTS  OF  UNFILTERED  AND  FILTERED 


IRSTD  OUTPUTS 


LMS  Filter 

alx  IIP 

X  W 

LMSI 

9.94 

14.9 

LMS2 

9.94 

19.3 

LMS3 

9.94 

29.3 

LMS4 

9.94 

39.3 

LMS5 

9.94 

49.3 

L\IS6 

9.94 

59.3 

LMS7 

9.94 

69.3 

LMS8 

9.94 

79.2 

L\1S9 

9.94 

89. 2 

LMSK) 

9.94 

99.2 

LMSll 

9.94 

109.1 

LMS12 

9.94 

119.0 

LMS13 

9.94 

128.9 

LMS14 

9.94 

138.8 

LMS15 

9.94 

148,6 

C.  FILTER  OUTPUTS 

Various  selections  of  the  15  trial  filters  were  used  on  BKDl,  BKDII,  and  BKDIIl. 
TGTl  through  TGTVI  were  embedded  into  each  background,  and  the  resultant  output 
signals  were  filtered.  The  S  N  ratios  of  both  the  filtered  and  unfiltered  outputs  were  then 
determined,  and  are  listed  in  Table  4  on  page  45,  Table  5  on  page  45,  and  Table  6  on 
page  45. 
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Table  4.  S/N  RATIOS  FOR  BKD  I 


TGT 

LMS  Filter 

ra'> 

LMS8 

LMS9 

LMSIO 

LMSll 

LMS  12 

LMS  13 

LMS14 

1 

0.932 

29.0 

29.7 

29.1 

27.9 

26.4 

24.8 

23.2 

11 

0.936 

29.0 

29.8 

29.2 

28.0 

26.5 

24.9 

23.3 

III 

0.934 

28. 3 

29.2 

28.8 

27.7 

26.4 

24.9 

23.3 

IV 

WKIHi 

23.5 

25.2 

25.8 

25.6 

25.0 

24.2 

23.1 

0.936 

19.1 

21.0 

22.1 

22.7 

22.9 

22.8 

22.4 

VI 

0.959 

18.7 

20.4 

21.4 

22.0 

22.4 

22.5 

22.4 

Table  5.  S/N  RATIOS  FOR  BKD  II 
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These  tables  show  that  L.MS9  through  LMS13  yield  the  highest  S  N  ratios  for  the 
simulated  backgrounds,  and  that  L.MS3  through  L.MS5  yield  the  highest  S  N  ratios  for 
the  real  background.  This  diflerence  between  backgrounds  is  expected,  because  the  real 
background  was  sampled  at  a  rate  three  times  slower  than  the  simulated  background. 
This  increased  sampling  time  corresponds  directly  to  a  reduction  in  the  half-amplitude 
width  to  sampling  time  ratio  ( -p )  by  a  factor  of  three.  Therefore,  taking  the  case  of 
TGTl  as  an  e.xample,  the  change  in  the  optimal  filter  from  LMS9  for  BKDl  and  BKDII 
to  L.MS3  for  BKDIII  is  expected. 

Additionally,  as  discussed  by  Niizberg  et  al.  (Ref.  3),  an  increase  in  the  sampling  time 
of  a  digital  system  leads  to  a  sampling-effect  loss  of  signal  shape  information  in  the 
system  output.  This  loss  of  signal  shape  information  can  be  considered  as  a  source  of 
high  frequency  noise,  which  will  be  amplified,  along  with  the  target  signal,  in  the  filtering 
process.  Therefore,  the  increased  sampling  time  of  BKDIII  accounts  for  the  remarkable 
difference  between  the  S  N  ratios  of  the  filtered  outputs  of  the  real  and  simulated  back¬ 
grounds. 

For  purposes  of  comparison,  the  reduction  of  S  N  ratio  with  increased  sampling 
times  for  the  LMS  filter  was  studied  by  Nitzberg  ei  al.  in  a  white  noise  environment  [Ref 
3J.  This  study  showed  a  decrese  in  S  N  ratio  by  a  factor  of  approximately  three,  as  the 
sampling  time  was  increased  by  a  factor  of  three.  The  results  of  this  study  compare  well 
with  the  decrease  in  S  N  ratio  by  a  factor  of  approximately  five  in  the  present  case. 
This  comparison  can  be  made,  because  the  study  by  Nitzberg  ei  al.  considered  only  a 
white  noise  background,  whereas  the  present  case  considers  a  background  consisting  of 
white  noise  as  well  as  background  clutter,  both  of  which  are  known  to  be  corrupted  b> 
quantization  error  due  to  digitization.  These  additional  noise  sources  are  amplified, 
along  with  the  white  noise,  as  the  sampling  time  is  increased  from  0.10472  mrad  sample 
to  0.31416  mrad  sample,  thereby  causing  additional  degradation  of  the  S/N  ratio.  Ad¬ 
ditionally.  since  the  IRSTD  DCU  samples  data  at  0.10472  mrad  sample,  vice  0.31416 
mrad  sample,  sampling  time  degradation  is  not  expected  from  IRSTD  outputs  digitized 
by  the  system's  DCU. 

To  illustrate  the  effectiveness  of  these  filters,  various  filtered  and  unfiltered  outputs 
are  illustrated  in  Figure  26  on  page  47  through  Figure  35  on  page  57  of  this  Chapter, 
as  well  as  in  Figure  36  on  page  61  through  Figure  61  on  page  78in  Appendix  .A  and 
Appendix  B.  The  significance  of  these  figures  will  be  discussed  in  the  following  sections. 
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1.  Unfiltered  Outputs 

Figure  26  on  page  47,  Figure  27  on  page  48,  and  Figure  28  on  page  48  are 
plots  of  the  unfiltered  outputs  of  BKDI,  BKDII,  and  BKDIII  with  TGTI,  TGTI,  and 
TGTVl  embedded,  respectively.  The  targets  are  located  at  the  center  of  the  unfiltered 
outputs,  and  have  a  S  N  ratio  of  approximately  one  in  BKDI  and  BKDll,  and  a  S  N 
ratio  of  approximately  two  in  BKDlll.  These  plots  are  shown  to  illustrate  the  fact  that 
extraction  of  the  targets  from  the  background  by  simple  means,  such  as  thresholding, 
would  not  be  very  effective. 

I  o  1 


in 


AZIMUTH  (MRRD) 


Figure  26.  Unfiltered  Output;  BKDI,  TGTI;  TGTI  embedded  in  BKDI; 
S  N  =  0.932. 
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in 


0.0  100. 0  200.0  300.0  400.0 

_ nZIMUTH  (HRRD) _ 

Figure  27.  Unfiltered  Output;  BKDII,  TGTl:  TGTI  embedded  in  BKDII; 
S,N  =  0.909. 


0.00  93.75  187.50  281.25  375.00 


_ RZIMUTH  (MRRD) _ 

Figure  28.  Unfiltered  Output;  BKDIII,  TCTVI:  TGTVI  embedded  in  BKDIII; 
S  N=1.95. 
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2.  Filtered  Outputs:  Target  Amplitude  versus  Noise  Amplitude 

Figure  29  on  page  51  through  Figure  32  on  page  54  are  plots  of  TGTl  embed¬ 
ded  in  BKDI  and  filtered  by  LMS  filters  of  increasing  width;  L.MS6,  LMS9,  LMS12,  and 
L\1S15  respectively.  It  can  be  seen,  by  referring  to  the  upper  plot  of  each  figure,  that 
the  amplitude  of  the  target  signal  increases  with  filter  width,  as  expected  from  the 
thresholding  properties  of  the  filters.  However,  the  S/N  ratio  reaches  a  maximum  when 
the  filter  s  half-amplitude  width  most  nearly  matches  the  half-amplitude  width  of  the 
target  signal  { Figure  30  on  page  52). 

.Additionally,  as  the  filter  widens  more  and  more  noise  is  also  allowed  to  pass, 
thereby  reducing  the  S  N  ratio.  This  increase  in  noise  with  filter  width  is  displayed  in 
the  lower  plot  of  each  figure,  which  has  been  normalized  to  an  amplitude  of  1.0  to  give 
a  better  illustration  of  the  noise.  It  can  be  seen  that  the  amplitudes  of  individual  noise 
peaks  begin  to  increase  rapidly  as  the  filter  width  increases  past  L.V1S9  (Figure  30  on 
page  52).  Specifically,  as  the  filter  widths  increase  to  LMS  12  and  LMS  15  (Figure  31 
on  page  53  and  Figure  32  on  page  54),  sharp  peaks  begin  emerging  from  the  background 
noise.  These  peaks  are  undesirable,  in  that  they  could  be  interpreted  as  false  targets  if 
the  post-filter  thresholding  were  set  too  low. 

From  studying  these  figures,  it  can  be  concluded  that  LMS9  results  in  the  opti¬ 
mal  output  for  BKDI.  This  filter  not  only  amplifies  TGTl  to  a  level  appropriate  for 
detection  and  acquisition  processing,  but  passes  a  minimal  anunount  of  noise,  compared 
to  filters  of  greater  width. 

Figure  36  on  page  61  through  Figure  39  on  page  64,  in  .Appendix  ,A.  arc  plots 
of  TGTl  embedded  in  BKDll  and  filtered  by  L.MS6,  LMS9,  L.MS12,  and  LMS15.  Fig¬ 
ure  40  on  page  65  through  Figure  43  on  page  6S,  in  Appendix  A,  are  plots  of  TGI  1 
embedded  in  BKDI  11  and  filtered  by  L.MSl,  LMS3,  L.MS5,  and  LMS7.  Using  similar 
analysis  on  these  figures,  it  can  be  concluded  that  LMS9  results  in  the  optimal  output 
for  TGTl  in  BKDII,  and  L.MS3  results  in  the  optimal  output  for  TGTl  in  BKDllI. 
However,  before  it  is  possible  to  consider  LMS3  and  LMS9  optimal  filters  for  other 
targets  in  these  backgrounds,  it  is  necessary  to  look  at  the  ability  of  these  filters  to  am¬ 
plify  TGTl  I  through  TGTVL 

Referring  back  to  Table  4  on  page  45,  Table  5  on  page  45,  and  Table  6  on  page 
45,  it  can  be  seen  that  TGTl  through  TGTHI  have  optimal  S/N  ratios  for  LMS3  and 
L.MS9.  whereas  TGTIV  through  TGTVl  have  optimal  S'N  ratios  for  LMS  filters  of 
greater  width.  However,  the  S  N  ratios  of  TGTIV  through  TGTVl  for  LMS3  and 
LMS9  arc  very  close  in  magnitude  to  the  maximum  S'N  ratios  for  these  targets,  within 
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16“o  in  the  worst  case.  Therefore,  it  can  be  assumed  that  LMS3  and  LMS9  will  amplify 
all  six  targets  to  levels  appropriate  for  post-filter  processing. 

To  illustrate  that  this  assumption  is  valid,  plots  of  all  target  and  background 
combinations,  filtered  with  these  proposed  optimal  filters,  are  included  as  Figure  44  on 
page  69  through  Figure  61  on  page  78  in  Appendix  B.  As  these  plots  illustrate,  LMS3 
and  LMS9  efi'ectively  remove  the  DC  ofl'set  and  linear  noise  from  the  background  clut¬ 
ter.  and  amplify  the  target  signals  for  all  background  and  target  combinations.  These 
resultant  outputs  can  then  be  effectively  used  as  inputs  for  target  acquisition,  desig¬ 
nation,  and  tracking  routines,  for  all  target  and  background  combinations  similar  to 
those  used  in  this  study. 
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Filtered  Output;  BKDI,  TGTI,  LMS6:  Upper  plot,  direct  filter  out¬ 
put;  lower  plot,  normalized  to  peak  amplitude  of  1.0;  S  'N'  =  26.7. 
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Figure  30.  Filtered  Output;  BKDI,  TGTI,  LMS9:  Upper  plot,  direct  filter  out¬ 
put;  lower  plot,  normalized  to  peak  amplitude  of  1.0;  S  N=  29.7. 


AMPLITUDE 

-G.5  0.0  0.5 


0.0  100.0  200.0  300.0  400.0 


AZIMUTH  (MRAD) 


o 


0.0  100.0  200.0  300.0  400.0 

AZIMUTH  (MRAD) 


Figure  31.  Filtered  Output;  BKDI,  TGTI,  LMS12:  Upper  plot,  direct  filter  out¬ 
put;  lower  plot,  normalized  to  peak  amplitude  of  1.0;  S  N  =  26.-4. 
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Figure  32.  Filtered  Output;  BKDI,  TGTI,  LMS15:  Upper  plot,  direct  filter  out¬ 


put;  lower  plot,  normalized  to  peak  amplitude  of  1.0;  S  N  =  21.6. 


3.  Filtered  Outputs:  Multiple  Targets 

Figure  33  is  a  plot  of  the  filtered  output  of  TGTI,  TGTIV,  and  TGTV  embed¬ 
ded  in  BKDI.  and  filtered  with  LMS9.  Figure  34  on  page  56  is  a  plot  of  the  filtered 
output  of  TGTI,  TGTIV,  and  TGTVI  embedded  in  BKDIII,  and  filtered  with  LMS3. 
As  these  plots  illustrate,  LMS3  and  L.VIS9  provide  efiective  filtering  of  multiple  targets 
of  vary  ing  size,  which  is  essential  if  they  are  to  be  used  in  a  real  environment. 

4.  Filtered  Outputs:  Real  Background  and  Target 

Figure  35  on  page  57  is  a  plot  of  real  IRSTD  output  data  obtained  during  sys¬ 
tem  calibration  [Ref  2].  The  upper  plot  shows  the  unfiltered  output  as  the  system  was 
scanned  across  the  roof  of  Spanagel  Hall  at  NFS  Monterey.  The  lower  plot  shows  the 
effect  of  L.V1S3  on  the  raw  system  output.  As  illustrated,  the  background  contains  many 
clutter  sources,  in  addition  to  sky.  and  LMS3  proves  to  be  extremly  effective  in  removing 
them  all.  as  well  as  filtering  out  the  whip  antenna  target  (TGTVI).  This  filtered  output 
clearly  illustrates  the  potential  of  the  LMS  filter  as  an  effective  signal  processing  filter 
for  the  IRSTD  system. 
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Figure  33.  Filtered  Output;  Multiple  Targets:  TGTI,  TGTIV,  and  TGTV  em¬ 
bedded  in  BKDI  and  filtered  with  LMS9. 
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Figure  34.  Filtered  Output;  Multiple  Targets:  TGTI,  TGTIV,  and  TGTVl  em¬ 
bedded  in  BKDIII  and  filtered  with  LMS3. 
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Figure  35.  Real  Output:  Upper  plot,  untiltered  output;  lower  plot,  filtered  with 
LMS3. 
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V  II.  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  CONCLUSIONS 

From  the  results  of  Chapter  VI,  it  can  be  assumed  that  LMS9  and  LMS3  are  the 
optimal  LMS  filters  for  extracting  point  and  near-point  targets  from  the  simulated  and 
real  I R  backgrounds  respectively.  This  diflerence  in  optimal  filters  for  the  simulated  and 
real  backgrounds  is  due  to  the  fact  that  the  simulated  background  was  sampled  at 
0.10472  mrad  sample,  whereas  the  real  background  was  sampled  at  0.31416 
mrad  sample.  Therefore,  these  optimal  filters  were  generated  from  a  raised  cosine  func¬ 
tion  with  a  half-amplitude  width  of  0.9  mrad,  sampled  at  0.1  mrad  sample  for  LMS9  and 
0.3  mrad  sample  for  LMS3.  This  result  corresponds  with  prefiltering  studies  on  the 
half-amplitude  width  of  the  target  signals,  and  is  an  indication  that  the  raised  cosine 
target  model  is  appropriate  for  these  targets. 

It  should  also  be  noted  that  the  results  of  this  investigation  can  be  used  with 
standard  inputs,  indepeiident  of  the  local  background,  whereas  other  approaches  [Refs. 
9,  13]  require  a  priori  knowledge  of  the  background  structure. 

Since  the  LMS  filter  is  ver}'  simple  and  takes  minimal  processing  time,  it  is  ideally 
suited  for  implementation  as  an  initial  filter  for  the  NFS  IRSTD  system.  Additionally, 
since  the  output  of  the  IRSTD  system  consists  of  digitized  detector  outputs  from  an 
array  of  90  detectors,  it  is  natural  to  consider  extending  this  filtering  technique  to  two- 
dimensions.  However,  since  the  output  signal  is  digitized  at  0.10472  mrad  per  sample  in 
the  azimuth  direction,  and  the  detector  dimensions  are  2.0  mrad  in  the  elevation  direc¬ 
tion.  the  resolution  of  the  two-dimensional  output  signal  can  be  considered  to  be  0.10472 
mrad  by  2.0  mrad  per  data  point.  Considering  this  order  of  magnitude  difference  in  the 
resolution  of  the  two  directions  of  the  two-dimensional  signal,  and  the  fact  that  the  point 
and  near-point  targets  studied  have  half-amplitude  widths  of  approximately  0.9  mrad,  it 
can  be  seen  that  filtering  in  the  elevation  direction  would  necessitate  the  use  of  an  LMS 
filter  generated  from  a  model  target  function  with  a  half-amplitude  width  to  sampling 
time  ratio  of  one  (-^  =  1).  Therefore,  a  suitable  two-dimensional  LMS  filter  model  tar- 
get  function  would  have  to  be  developed  from  a  two-dimensional  Gaussian  distribution 
or  a  raised  Bessel  function  defined  over  its  first  period  of  oscillation.  In  addition,  either 

model  target  function  would  have  to  be  constructed  in  such  a  manner  that  -p  =  9  in  the 
^  t 

azimuth  direction  and  -p  =  1  in  the  elevation  direction. 
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Although  developing  and  implementing  such  a  filter  is  feasible,  the  processing  time 
required  for  performing  90  one-dimensional  filtering  operations  over  N  elements,  is 
comparable  with  the  processing  time  required  for  performing  one  two-dimensional  fil¬ 
tering  operation  over  90  by  N  elements  (2.14  vice  2.03  seconds  respectively).  Therefore, 
it  is  advisable  to  implement  this  filtering  technique  in  its  simplest,  or  one-dimensional, 
form. 

B.  RECOMMENDATIONS 

It  is  recommended  that  the  following  studies  be  conducted  on  the  NFS  IRSTD  sys¬ 
tem: 

•  Incorporate  LMS9  as  an  initial  signal  processing  filter  on  the  output  of  the  lead 
detector  array  of  the  IRSTD  system. 

•  Follow  LMS9  with  an  automatic  target  acquisition  routine. 

•  Study  the  outputs  of  this  automatic  target  acquisition  routine  to  optimize  target 
false  alarm  rates. 

•  Study  the  outputs  of  various  hypothesis  testing  techniques  (Ref  14  ],  to  further 
discriminate  between  true  targets  (missiles  and  aircraft),  and  false  targets  (birds, 
high  frequency  noise,  background  discontinuities,  etc.). 

•  Develop  a  filtering  technique  for  extracting  extended  targets  (t.,>  1.5  mrad)  from 
background  clutter. 

•  Incorporate  this  extended  target  filter  into  the  IRSTD  system,  to  filter  the  outputs 
from  the  lag  detector  array.  This  would  give  the  system  two  dilTerent  outputs  for 
each  scan. 

•  Study  the  potential  for  imaging  the  extended  targets,  with  the  intention  of  devel¬ 
oping  ship-type  recognition  techniques. 

It  is  believed  that  much  groundwork  for  further  studies  with  the  NFS  IRSTD  system 
has  been  laid  by  this  thesis.  Follow-on  studies  of  the  types  listed  above  are  both  feasible 
and  important,  as  the  U  S.  Na\T  has  great  need  for  an  effective  passive  IR  system  ca¬ 
pable  of  detecting  and  tracking  air  and  surface  targets. 


APPENDIX  A.  FILTERED  OUTPUTS:  TARGET  AMPLITUDE  VERSUS 

NOISE  AMPLITUDE 


Figure  36  on  page  61  through  Figure  43  on  page  68  are  plots  of  filtered  outputs  of 
TGTI  embedded  in  BKDII  and  BKDllI.  BKDll  was  filtered  by  LMS6.  LMS9,  LMS12, 
and  LMS15.  and  BKDIII  was  filtered  by  LMSl,  LMS3,  LMS5,  and  LMS7.  It  can  be 
seen,  by  referring  to  the  upper  plot  of  each  figure,  that  the  amplitude  of  the  target  signal 
increases  with  filter  width,  as  expected  from  the  thresholding  properties  of  the  filters. 
However,  the  S  N  ratio  reaches  a  maximum  when  the  filter's  half-amplitude  width  most 
nearly  matches  the  half-amplitude  width  of  the  target  signal  (Figure  37  on  page  62  and 
Figure  41  on  page  66). 

Additionally,  as  the  filler  widens  more  and  more  noise  is  also  allowed  to  pass, 
thereby  reducing  the  S  N  ratio.  This  increase  in  noise  with  filter  width  is  displayed  in 
the  lower  plot  of  each  figure,  which  has  been  normalized  to  an  amplitude  of  1 .0  to  give 
a  better  illustration  of  the  noise.  It  is  clear,  from  these  figures,  that  the  noise  amplitude 
increases  rapidly  with  filter  width,  and  contains  sharp  peaks  which  could  be  interpreted 
as  false  targets  if  post-filter  thresholding  were  set  too  low. 

From  studying  these  figures,  it  can  be  concluded  that  L.V1S9  results  in  the  optimal 
output  for  BKDII.  and  LMS3  results  in  the  optimal  output  for  BKDIII.  These  filters 
not  only  amplify  TGTI  to  levels  appropriate  for  detection  and  acquisition  processing, 
but  pass  a  minimal  ammount  of  noise,  compared  to  filters  of  greater  width. 
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Filtered  Output:  BKDII,  TGTI,  LMS6:  Upper  plot,  direct  filter  out¬ 
put;  lower  plot,  normalized  to  peak  amplitude  of  I.O;  S;N  =  28.9. 


Figure  37. 


Filtered  Output;  BKDII,  TGTI,  LMS9:  Upper  plot,  direct  filter  out¬ 
put;  lower  plot,  normalized  to  peak  amplitude  of  1.0;  S/N*  32.2. 
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Figure  38.  Filtered  Output;  BKDII,  TGTI,  LMS12:  Upper  plot,  direct  filter 
output;  lower  plot,  normalized  to  peak  amplitude  of  1.0;  S  N  =  28.0. 
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Figure  39.  Filtered  Output;  BKDII,  TGTI,  LMSI5:  Upper  plot,  direct  filter 
output;  lower  plot,  normalized  to  peak  amplitude  of  1.0;  S,'N  =  21.6. 
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Figure  40.  Filtered  Output;  BKDIIU  TGTI,  LMSl:  Upper  plot,  direct  filter 
output;  lower  plot,  normalized  to  peak  amplitude  of  1.0;  S/N*  8.1. 


65 


AMPLITUDE 

“0.5  0.0  0.5 


1 

iT 

1 

D- 

1. 

93.75  187.50  281.25 

AZIMUTH  (MRAD) 


375.00 


93.75  187.50  281.25 

AZIMUTH  (MRAD) 


375.00 


66 


AMPLITUDE 

.50-0.25  0.00  0.25  0.50 


67 


Figure  43. 


Filtered  Output;  BKDIII,  TGTI,  LMS7:  Upper  plot,  direct  filter 
output;  lower  plot,  normalized  to  peak  amplitude  of  1.0;  S,’N=  7.6. 
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APPENDIX  B.  FILTERED  OUTPUTS:  OPTIMAL  LMS  FILTERS 


Figure  44  through  Figure  61  on  page  78  are  plots  of  the  filtered  outputs  of  all 
background  and  target  combinations  filtered  with  their  proposed  optimal  filters,  LMS3 
and  LMS9.  As  these  plots  illustrate,  these  filters  effectively  remove  the  DC  offset  and 
linear  noise  from  the  background  clutter,  and  amplify  the  target  signals  for  all  back¬ 
ground  and  target  combinations.  These  resultant  outputs,  can  be  effectively  used  as  in¬ 
puts  for  target  acquisition,  designation,  and  tracking  routines. 


Figure  44.  Filtered  Output;  BKDI,  TGTI,  LMS9:  TGTl  embedded  in  BKDl  and 

filtered  by  LMS9;  S  N  =  29.7. 


69 


Figure  45.  Filtered  Output;  BKDI,  TGTII,  LMS9:  TGTII  embedded  in  BKDI 
and  filtered  by  LMS9;  S;N=  29.8. 
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Figure  46.  Filtered  Output;  BKDI,  TGTIII,  LMS9:  TGTIII  embedded  in  BKDI 
and  filtered  by  LMS9;  S,'N  =  29.2. 
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Figure  47.  Filtered  Output;  BKDl,  TGTIV,  LMS9:  TGTIV  embedded  in  BKDI 
and  filtered  by  LMS9;  S/N=  25.2. 


Figure  48.  Filtered  Output;  BKDI,  TGTV,  LMS9:  TGTV  embedded  in  BKDI 
and  filtered  by  LMS9;  S/N=  21.0. 
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Figure  49.  Filtered  Output;  BKDI,  TGWI,  LMS9:  TGTVI  embedded  in  BKDI 
and  filtered  by  LMS9;  S/X  =  20.4. 


Figure  50.  Filtered  Output;  BKDII,  TGTI,  LMS9;  TGTI  embedded  in  BKDII 
and  filtered  by  LMS9;  S/N=  32.2. 
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Filtered  Output;  BKDII,  TGTII,  LMS9:  TGTII  embedded  in  BKDII 
and  filtered  by  LMS9;  S.X  =  32.3. 
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Filtered  Output;  BKDII,  TGTIII,  LMS9:  TGTII  I  embedded  in 
BKDII  and  filtered  by  LMS9;  S,X=3I.7. 
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Figure  53.  Filtered  Output;  BKDII,  TGTIV,  LMS9:  TGTIV  embedded  in 
BKDII  and  filtered  by  LMS9;  S  N=27.3. 
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Figure  55.  Filtered  Output;  BKDII,  TGTVI,  LMS9:  TGTVI  embedded  in 
BKDII  and  filtered  by  LMS9;  S;N=22.1. 
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Figure  56.  Filtered  Output;  BKDIII,  TGTI,  LMS3:  TGTI  embedded  in  BKDIII 
and  filtered  by  LMS3;  S.N=  12.7. 
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57.  Filtered  Output;  BKDIIl,  TGTII,  LMS3:  TGTII  embedded  in 
BKDIII  and  filtered  by  LMS3;  S;N=  12.7. 
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Figure  58.  Filtered  Output;  BKDIII,  TGTIII,  LMS3;  TGTIII  embedded  in 
BKDIII  and  filtered  by  LMS3;  S  N=  12.5. 
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Figure  59.  Filtered  Output:  BKDIII,  TGTIV,  LMS3:  TGTIV  embedded 
BKDIII  and  filtered  by  LMS3;  S  N  =  11.1. 
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Figure  60.  Filtered  Output;  BKDIII,  TGTV,  LMS3:  TGTV  embedded  in 
BKDIII  and  filtered  by  LMS3;  S/N  =  9.4. 
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Figure  61.  Filtered  Output;  BKDIII,  TGTVI,  LMS3:  TGTVI  embedded  in 
BKDIII  and  filtered  by  LMS3;  S.N  =  9.1. 
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APPENDIX  C.  BACKGROUND  MODEL  FORTRAN  PROGRAM 


This  appendix  is  included  to  present  the  FORTRAN  code  used  for  generating  sim¬ 
ulated  background  distributions.  All  comments  and  discussion  of  the  parameters  and 
functions  of  the  program  are  included  as  comment  lines  in  the  program  code. 

This  program  was  designed  to  be  run  on  the  NTS  IBM  3033/4381  Mainframe 
Computer  Network,  using  VS  FORTR*\N  on  the  systems  Virtual  Machine.  The  output 
of  this  program  contains  the  simulated  IR  background  amplitude  distribution  in  a  20 
by  4000  element  matrix. 


C  PROGRAM  BKDl 

DIMENSIO.M  U1(4000,20),U2(201),U3(4Q00),U4(101) 
DOUBLE  PRECISION  DSEED 
CHARACTER*80  T1$,XL$ 

CHARACTER*40  CL$ , SI$ 

CALL  SHERPA  ( ' BKDPL0T2 ' , ' A' , 3) 

C 

C  CO.NSTANTS  AND  VARIABLES 

C 

M=4000 

M2=201 

M3=101 

N=20 

IT=1 

DSEED=414813567. DO 
CLEN=1667. 

SIGMA=0. 15 

CLS=' CORRELATION  LENGTH  =  166.7  MRAD?' 

SIS=' SIGMA  =  0.  1$' 

IZ=INT( DSEED) 

ST=SQRT(2.  ) 

Al=-1. /CLEN 
Z=EXP(A1) 

Z2=Z/(1.+(Z**ST)) 

Bl=l.  -(Z*Z) 

B2=SQRT(B1) 

B3=l.  -(2.*Z*Z)/(1.+(Z**ST)) 

B4=SQRT(B3) 

DO  2  1=1, N 
DO  i  J=1,M 
UKJ,I)=0. 

1  CONTINUE 

2  CONTINUE 
C 

C  SET  FIRST  ELEMENT  EQUAL  TO  ZERO 
C 
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1000  U1(1,1)=0. 

GENERATE  RANDOM  NUMBERS  FOR  FIRST  ROW 

CALL  GGNML(DSEED,M,U3) 

GENERATE  ELEMENTS  OF  FIRST  ROW 

DO  20  1=2, M 
11=1-1 
RN=U3(I1) 

UK  1 , 1  )=(  SIGMA*B2*RN)+(  Z*U1(  1 1 , 1)  ) 

20  CONTINUE 

GENERATE  REST  OF  BACKGROUND  (ROW  BY  ROW) 

DO  50  1=2, N 

GENERATE  RANDOM  NUMBERS  (ROW  BY  ROW) 

CALL  GGNML(DSEED,M,U3) 

DO  40  J=1,M 
11=1-1 
J1=J-1 
RN=U3(J) 

GENERATE  ELEMENTS  OF  FIRST  COLUMN 
IF  (J.EQ.  1)  U1(J,I)=(SIGMA*B2*RN)+(Z*U1(J,I1)) 
GENERATE  REST  OF  ELEMENTS 

U1(J,I)=(SIGMA*B4*RN)+(Z2*(U1(J,I1)+U1(J1,I))) 
40  CONTINUE 
50  CONTINUE 

PLOT  3D  BACKGROUND  DISTRIBUTION 

CALL  PL0T32  (U1,M,N) 

DETERMINE  AMPLITUDE  DISTRIBUTION 

CALL  DIST  (U1,U2,U4,M,M2,M3,N,IT,CL$,SI$,IZ) 

IF  (IT.  LE.  1)  THEN 
IZ=INT(DSEED) 

GO  TO  1000 
END  IF 
CALL  DONEPL 
STOP 
END 


SUBROUTINE  DIST(U1 ,U2 ,U3 ,M,M2,M3,N, IT,CL$ ,SI$ , IZ) 
DIMENSION  U1(M,N),U2(M2),U3(M3) 

CHARACTER*40  CL$,SI$ 

DO  5  J=1,M2 
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U2(J)=0. 

5  CONTINUE 
DO  20  1=1, N 
DO  10  J=1,M 

II=INT(U1(J,I)*100.  )+101. 
U2(II)=U2(II)+1. 

10  CONTINUE 
20  CONTINUE 
YMAX=0. 

DO  40  1=1, M2 

IF  (U2(I).GT.  YMAX)  YMAX=U2(I) 
40  CONTINUE 

DO  50  I=51.M2-50 

J=M2-50+l-I 

U3(J)=U2(I) 

50  CONTINUE 
DO  60  1=1, M2 
U3(I)=U3(I)/YMAX 
60  CONTINUE 

PLOT  AMPLITUDE  DISTRIBUTION 

CALL  PLOTIT(U3,M3) 

IT=IT+1 

RETURN 

END 


SUBROUTINE  PLOTIT  (Y,N) 

DIMENSION  Y(N),X(101) 

REAL  X,Y 
INTEGER  IFIG 
CALL  PRTPLT(76,2) 

CALL  TEK618 

CALL  PAGE  (11.  ,8.5) 

CALL  NOBRDR 

CALL  HVROT( 'MOVIE' ) 

CALL  AREA2D  (4. 8,2. 6) 

CALL  HEIGHT  (0. 15) 

CALL  XNAME  ( ' AMPLITUDE$ ' , 100) 

CALL  YNAME  ('OCCURANCES  (NORMALIZED) $’ ,100) 

55  CALL  GRAF  ( -,  5 , ’ SCALE ' , .  5 ,0.  , ’ SCALE' , 1.  ) 

DO  60  1=1, N 

X(I)=. 01*FLOAT(I)-.51 

60  CONTINUE 

65  CALL  CURVE  (X,Y,N,0) 

CALL  ENDPL(O) 

RETURN 

END 


SUBROUTINE  PL0T32(U1 ,M,N) 
DIMENSION  U1(M,N) 
AN=REAL(N)*.  1 
AM=REAL(M)*.  1 
AN 1= AN/ 2. 
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AMl=AM/4. 

CALL  PRTPLT(76,2) 

CALL  TEK618 

CALL  PAGE  (11.  ,8.  5) 

CALL  NOBRDR 

CALL  HWROT( 'MOVIE') 

CALL  AREA2D  (5.  ,3.  ) 

CALL  HEIGHT  (0. 15) 

CALL  X3NAME  ('AZIMUTH  (MRAD)$ ' , 100) 

CALL  Y3NAME  ('ELEVATION  (MRAD)$ ' , 100) 

CALL  Z3NAME  ( 'AMPLITUDE$' ,100) 

CALL  V0LM3D  (2. ,2.  ,1.  ) 

CALL  VUABS  (5.  ,-S.  ,3.  ) 

CALL  GRAF3D  ( 0.  ,AM1 , AM,0.  ,AN1,AN,-.  5,.  2,.  5) 
CALL  SURMAT  (U1,M,M,0,N,0) 

CALL  ENDPL(O) 

RETURN 

END 
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APPENDIX  D.  IRSTD  MODEL  FORTRAN  PROGRAM 


This  appendix  is  included  to  present  the  FORTRAN  code  used  for  generating  sim¬ 
ulated  IRSTD  detector  outputs.  All  comments  and  discussion  of  the  parameters  and 
functions  of  the  program  are  included  as  comment  lines  in  the  program  code. 

This  program  was  designed  to  be  run  on  the  NPS  IBM  3033  4381  Mainframe 
Computer  Network,  using  VS  FORTRAN  on  the  systems  Mutilple  Virtual  System. 
The  output  of  this  program  contains  the  simulated  IR  detector  output  in  a  4000  element 
array. 


//GRIBIIRS  JOB  (3222, 9999), ’ IR1D1MVS6' ,CLASS=G,REGION=3000K 
//*F0RMAT  PR,DDNA.ME=G0.  FT06F001, 

//*F0RMS=SEP1 

//  EXEC  F0RTVCLG,PARM.  F0RT='0PT(2)' 

//FORT. SYS IN  DD  * 

C 

C  IRST  MODEL  -  (4000,20)  OBJECT  PLANE 
C 

IMPLICIT  COMPLEX  (C) 

DI.ME.NSION  01(4000,20)  ,U2(4000), 03(4000) 

DIMENSION  06(4000) ,07(4000) ,08(4000) 

DIMENSION  01(4000,20) ,02(4000,20) 

DIMENSION  IWK(25150) ,RWK(25150) ,CWK(20) 

DOOBLE  PRECISION  DSEED 

DSEED=765742451.D0 

M=4000 

N=20 

IRAN=13455 

KRA.N=99947 

C 

C  GENERATE  BACKGROOND  DISTRIBUTION 

C 

CALL  BKGDl  (01 ,07 ,M,N,DSEED) 

DO  20  1=1, N 
DO  10  J=1,M 

C1(J,I)=CMPLX(01(J,I),0.  ) 

10  CONTINUE 
20  CONTINUE 
C 

C  GENERATE  SYSTEM  PSF 

C 

DO  40  1=1, N 
DO  30  J=1,M 

Al=REAL((J-(M/2))**2+(I-(N/2))**2)/40. 

IF  (Al.  GT.  100.  )  Al=100. 

A2=l.  /EXP(Al) 
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C2(J,I)=CMPLX(A2,0.  ) 
30  CONTINUE 
40  CONTINUE 


GENERATE  DETECTOR  ARRAY  DISTRIBUTION 

DO  50  J=1,M 
U2(J)=0. 

50  CONTINUE 
Jl=M/2 
J2=(M/2)+2 
DO  70  J=J1,J2 
U2(J)=1. 

70  CONTINUE 

GENERATE  IMAGE  PLANE  DISTRIBUTION 

CALL  FFT3D  (Cl ,M,N,M,N, 1 , 1 ,IWK,RWK,CWK) 
CALL  FFT3D  (C2 ,M,N,M,N, 1 , 1 , IWK,RWK,CWK) 
DO  120  1=1, N 
DO  110  J=1,M 
C1(J,I)=C1(J,I)*C2(J,I) 
no  CONTINUE 
120  CONTINUE 

CALL  FFT3D  (Cl ,M,N,M,N, 1 , -1 , IWK,RWK,CWK) 
DO  140  1=1, N 
DO  130  J=1,M 
U1(J,I)=CABS(C1(J,I)) 

130  CONTINUE 
140  CONTINUE 

GENERATE  IMAGE  PLANE  SLICES 

DO  180  J=1,M 
SUM 1=0. 

DO  170  Jl=l,20 
SUM1=SUM1+U1(J,J1) 

170  CO.NTINUE 
U6(J)=SUM1 
180  CONTINUE 

CALL  CONVOL  (U6,U2,U3,M) 

DO  230  J=1,M 
IRAN=IRAN*KRAN 

R=. 002*(REAL(IRAN)*2.  328306E-10) 
U3(J)=U3(J)+R 
230  CONTINUE 

AC  COUPLED  DETECTOR  OUTPUTS 

DO  260  J=1,M 
A1=REAL(J)/1.  3 
IF  (Al.  GT.  100.  )  GO  TO  250 
U7(J)=1.  /EXP(Al) 

GO  TO  260 
250  U7(J)=0. 

260  CONTINUE 
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CALL  CONVOL  (U3,U7,U8,M) 

WT1=0. 

DO  270  J=1,M 
U3(J)=U3(J)-U8(J) 

IF  (ABS(U3(J)).GT.WT1)  WT1=ABS(U3( J) ) 

270  CONTINUE 

DO  290  J=1,M 
U3(J)=.  35*(U3(J)/WT1) 

290  CONTINUE 

STORE  DETECTOR  OUTPUTS  IN  EXTERNAL  FILE 

DO  310  J=1,M,10 
WRITE(6,9000)(U3(J+K),K=0,9) 

9000  FORMAT  ( 10F7. 4) 

310  CONTINUE 
STOP 
END 
C 

SUBROUTINE  BKGD1(U1 ,U3 ,M,N,DSEED) 

DIMENSION  U1(M,N) ,U3(M) 

DOUBLE  PRECISION  DSEED 
DO  2  1=1, N 
DO  1  J=1,M 
U1(J,I)=0. 

1  CONTINUE 

2  CONTINUE 
CLEN=1623.  4 
SIGMA=0.  1 

C  CONSTANTS  AND  VARIABLES 
ST=SQRT(2. ) 

Al=-1.  /CLEN 
Z=EXP(A1) 

Z2=Z/(1.+(Z**ST)) 

Bl=l.  -(Z*Z) 

B2=SQRT(B1) 

B3=l. -(2.  *Z*Z)/( 1.  +(Z**ST)) 

B4=SQRT(B3) 

1000  U1(1,1)=0. 

CALL  GGNML( DSEED, M,U3) 

DO  20  1=2, M 
11=1-1 
RN=U3(I1) 

U1(I,1)=(SIGMA*B2*RN)+(Z*U1(I1,1)) 

20  CONTINUE 
DO  50  1=2, N 
CALL  GGNML( DSEED, M,U3) 

DO  40  J=1,M 
11=1-1 
J1=J-1 
RN=U3(J) 

IF  (J. EQ. 1)  U1(J,I)=(SIGMA*B2*RN)+(Z*U1(J,I1)) 
U1(J,I)=(SIGMA*B4*RN)+(Z2*(U1(J,I1)+U1(J1,I))) 
40  CONTINUE 
50  CONTINUE 
RETURN 
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END 

C 

SUBROUTINE  CONVOL(F,G,H,M) 
DIMENSION  F(M),G(M),H(M) 

DO  10  1=1, M 
SUM=0. 0 
DO  5  J=1,M 
K=M0D(I-J+M,M)+1 
SUM=SUM+F(J)*G(K) 

5  CONTINUE 
H(I)=SUM 
10  CONTINUE 
RETURN 
END 

/* 

// 
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APPENDIX  E.  FILTER  OUTPUT  FORTRAN  PROGRAM 

This  appendix  is  included  to  present  the  FORTRAN  code  used  for  filtering  simu¬ 
lated  and  real  background  distributions.  All  comments  and  discussion  of  the  parameters 
and  functions  of  the  program  are  included  as  comment  lines  in  the  program  code. 

This  program  was  designed  to  be  run  on  the  NPS  IBM  3033/4381  Mainframe 
Computer  Network,  using  VS  FORTRAN  on  the  systems  Virtual  Machine.  The  output 
of  this  program  contains  the  unfiltered  and  filtered  detector  outputs  in  5000  element  ar¬ 
rays.  and  the  results  of  S  N  ratio  calculations. 


C  PROGRAM  FILTl 

IMPLICIT  COMPLEX  (C) 

DIMENSION  U1(5000),U2(5000),U3(5000),XA(5000) 
DIMENSION  T(100),F(50) 

DIMENSION  C1(5000,1),C2(5000,1),C3(5000,1) 

DIMENSION  IVK( 30150), RWK( 30150), CWK(l) 

CHARACTER*80  TITLES, TIT$ 

CHARACTER*20  FILNS 

CALL  SHERPA  ( '0UTFILT2’ , 'A' ,3) 

C 

C  INITIALIZE  CONSTANTS 

C 

PI=3. 1415926535897932 

IT=1 

N2=100 

6000  DO  20  1=1,50 
F(I)=0. 

20  CONTINUE 
C 

C  SELECT  BACKGROUND 

C 

WRITEC*,*)' CHOOSE  BACKGROUND  (ENTER  BKD  NUMBER  1-5)’ 
READ( *)IBKD 

IF  (IBKD.EQ. 1)  FILN$='BKD1’ 

IF  (IBKD.EQ. 2)  FILN$='BKD2’ 

IF  (IBKD.EQ. 3)  FILN$='BKD3' 

IF  (IBKD.EQ.  4)  FILN$='BKD4’ 

IF  (IBKD.EQ.  5)  FILN$='BKD5’ 

N=4000 

IF  (IBKD.EQ. 3)  N=1250 
CALL  INPUT(U1,N, FILNS) 

C 

C  SELECT  TARGET 
C 

DO  10  1=1, N 
U2(I)=0. 
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10  CONTINUE 

WRITE(*,*)’ SELECT  NUMBER  OF  TGTS' 

READ(*,*)NTGT 
DO  35  I=1,NTGT 

WRITE(*,*)' CHOOSE  TARGETS  (ENTER  TGT  NUMBER  1-6)’ 
READ(*,*)ITGT 
IF  (IBKD.GE. 3)  GO  TO  25 
IF  (ITGT.EQ.  1)  FILN$='T1’ 

IF  (ITGT.EQ.  2)  FILN$='T2’ 

IF  (ITGT.EQ.  3)  FILN$='T3' 

IF  (ITGT.EQ.  4)  FILN$=’T4’ 

IF  (ITGT.EQ.  5)  FILN$='T5’ 

IF  (ITGT.EQ.  6)  FILN$='T6’ 

GO  TO  26 

25  IF  (ITGT.EQ.  1)  FILN$='T1A' 

IF  (ITGT.EQ.  2)  FILN$='T2A' 

IF  (ITGT.EQ. 3)  FILN$='T3A' 

IF  (ITGT.EQ.  4)  FILN$='T4A' 

IF  (ITGT.EQ. 5)  FILN$=’T5A’ 

IF  (ITGT.EQ.  6)  FILN$='T6A' 

26  CALL  INPUT(T,N2,FILN$) 

IF  (NTGT.  EQ.  1)  ITL=N/2 
IF  (NTGT.  NE.  1)  THEN 

WRITE(*,*)' ENTER  TARGET  LOCATION  (PIXELS  IN  INTEGER)' 

READ(*,*)ITL 

END  IF 

TL=0.  1*REAL(ITL) 

WRITE(*,*) 'ENTER  INTENSITY  FACTOR  OF  TARGET  (REAL)’ 

WRITE(*,*)'FOR  S/N=l:  0.15  BKDI  AND  BKDII;  0.09  BKDIIl’ 

READ(*,*)TFAC 

Kl=ITL-N2/2+l 

K2=K1+N2-1 

DO  30  J=K1,K2 

L=J-K1+1 

U2(J)=U2(J)+( TFAC*T( L) ) 

XA(J)=U2(J)+U1(J) 

IF  (XA(J).GT. 0.4095)  U2( J)=0.  4095-Ul( J) 

IF  (XA(J).LT. -0.4095)  U2( J)=-0.  4095-Ul( J) 

30  CONTINUE 
35  CONTINUE 
DO  40  J=1,N 

C1(J,1)=CMPLX(U1(J),0. ) 

40  CONTINUE 
DO  50  J=1,N 

C2(J,1)=CMPLX(U2(J),0.  ) 

50  CONTINUE 

SELECT  FILTER 

DO  60  J=1,N 
C3(J,1)=CMPLX(0.  ,0.  ) 

60  CONTINUE 

WRITE(*,*)' SELECT  FILTER  (ENTER  WIDTH-TO-SAMPLE  RATIO)’ 

READ(*,*)IFILT 

SUM=0. 

IFIL=(2*IFILT)+1 
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IFIL2=IFILT+1 
DO  70  I*1,IFIL 
II=I-IFIL2 

ARG=(PI*REAL( II))/REAL( IFILT) 
F(I)=COS(ARG)+l. 

SUM=  SUM+F(I) 

70  CONTINUE 

AVE=SUM/REAL( IFIL) 

AMAX=0. 

DO  80  1=1, IFIL 
F(I)=F(I)-AVE 

IF  (F(I).GT.AMAX)  AMAX=F(I) 

80  CONTINUE 

DO  90  1=1, IFIL 
F(I)=(F(I)/AMAX) 

90  CONTINUE 

DO  150  1=1, IFIL 

II=I-IFIL2 

K=N/2+II 

C3(K,1)=CMPLX(F(I),0.  ) 

150  CONTINUE 

FFT  OF  BKD,  TGT,  AND  FILTER 

CALL  FFT3D  (Cl ,N, 1 ,N, 1 , 1 , 1 ,IWK,RWK,CWK) 
CALL  FFT3D  (C2 ,N, 1 ,N, 1 , 1 , 1 , IWK,RWK,CWK) 
CALL  FFT3D  (C3,N, 1 ,N , 1 , 1 , 1, IWK,RWK,CWK) 

DATA  REDUCTION 

V=l.  /(12.*(2.**24)) 

VAR=2.  *V 
AS1=0. 

AN1S=0. 

AS2=0. 

AN2S=0. 

DO  160  1=1, N 
A1=CABS(C1(I,1)) 

A2=CABS(C2(I,1)) 

A3=CABS(C3(I,1)) 

AS1=AS1+A2 

AN1S=AN1S+(A1**2) 

AS2=AS2+(A2*A3) 

AN2S=AN2S+((A1*A3)**2) 

160  CONTINUE 

AN1=SQRT(AN1S) 

AN2=SQRT(AN2S) 

SN1=AS1/AN1 

SN2=AS2/AN2 

SH=0. 

DO  165  1=1, IFIL 
SH=SH+(F(I)**2) 

165  CONTINUE 
VARF=VAR*SH 
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C  EMBED  TARGET 
C 

KK=0 

DO  170  1=1, N 
IF  (U2(I).EQ.O.  )  THEN 
U1(I)=U1(I)+U2(I) 

GO  TO  170 
END  IF 

IF  (KK.EQ. 1)  GO  TO  166 
IF  (U2(I).NE.  0.  )  THEN 
KK=1 
BL=U1(I) 

END  IF 

166  U1(I)=U2(I)+BL 
170  CONTINUE 

DO  180  J=1,N 
C1(J,1)=CMPLX(U1(J).0.  ) 

180  CONTINUE 

FFT  OF  TGT  AND  BACKGROUND 

CALL  FFT3D  (C1,N,1,N,1,1,1,IWK,RWK,CWK) 

CONVOLVE  OUTPUT  WITH  FILTER 

DO  190  J=1,N 
C2(J,1)=C1(J,1)*C3(J,1) 

190  CONTINUE 

IFT  OF  FILTERED  OUTPUT 

CALL  FFT3D  (C2,N,1,N,1,1,-1,IWK,RWK,CWK) 

UNVRAP  OUTPUT 

DO  200  J=l,N/2 
K=J+N/2 

U3(J)=REAL(C2(K,1)) 

U3(K)=REAL(C2(J,1)) 

200  CONTINUE 

TITLE$=':  DETECTOR  OUTPUT  -  BKD  $’ 

WRITEC*,*)  'ENTER  1  IF  PLOT  OF  UNFILTERED  OUTPUT  DESIRED' 
READ(*,*)I0UT 
IF  (lOUT.  NE.  1)  GO  TO  210 

CALL  PLOTIT  (U1,XA,N, IT, TITLES, IBKD,ITGT,TL,SN1, VAR) 

IT=IT+1 

210  TIT$=' FILTERED  OUTPUT  $' 

CALL  PLOTI2  (U3,XA,N, IT, TITLE?, TIT?, IBKD,ITGT,IFILT,TL,SN2,VARF) 
IT=IT+1 

WRITEC*,*) 'ENTER  1  IF  ANOTHER  RUN  IS  DESIRED’ 

READ(*,*)IBREP 
IF  (IBREP.EQ,  1)  GO  TO  6000 
CALL  DONEPL 
STOP 
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END 


SUBROUTINE  INPUT(RE,N,FILN$) 
DIMENSION  RE(N) 

CHARACTER*20  FILN$ 

CLOSE  (2) 

OPEN( 2 ,FILE=FILN$ , STATUS=’ OLD’ ) 

DO  1  J  =  1,N,10 

READ( 2 , ' ( 10F7.  4) ' ) ( RE( J+K) ,K»0 , 9 ) 
1  CONTINUE 
RETURN 
END 


SUBROUTINE  PLOTIT  (Y,X,N,IFIG,TITLE$,IB,IT,TL,SN,VAR) 
DIMENSION  Y(N),X(N) 

CHARACTER*80  TITLE $ 

REAL  X,Y,STD 
INTEGER  IFIG 

IF  (IB.LE.  2)  AN=0.  1*REAL(N) 

IF  (IB.GT.  2)  AN=0.  3*REAL(N) 

ANS=AN/4. 

CALL  PRTPLT(76,2) 

CALL  TEK618 

CALL  PAGE  (11. ,8.5) 

CALL  NOBRDR 

CALL  HVROTC 'MOVIE') 

CALL  AREA2D  (5. ,2.5) 

CALL  HEIGHT  (0. 15) 

CALL  Y.VA.ME  ( 'AMPLITUDE?' ,100) 

CALL  XNAME  ('AZIMUTH  (MRAD)?’ ,100) 

CALL  MESSAG  ('FIGURE  $',100,1. ,5.6) 

CALL  INTNO  ( IFIG, 'ABUT' , 'ABUT' ) 

CALL  MESSAG  (TITLE$ , 100, 'ABUT' , 'ABUT' ) 

CALL  INTNO  (IB, 'ABUT' , 'ABUT' ) 

CALL  MESSAG  ('TARGET  NUMBER  $’,100,2.48,5.4) 

CALL  INTNO  ( IT, 'ABUT' , 'ABUT' ) 

CALL  MESSAG  ('TARGET  LOCATION  =  $’,100,2.48,5.2) 

CALL  REALNO  (TL, 1 , 'ABUT' , 'ABUT' ) 

CALL  MESSAG  ('  MRAD  $', 100, 'ABUT' ,’ ABUT' ) 

CALL  MESSAG  ('S/N  =  $’,100,2.48,5.0) 

CALL  REALNO  (SN,5, 'ABUT' , 'ABUT' ) 

CALL  MESSAG  ('VARIANCE  =  S’ ,100,2. 48,4. 8) 

CALL  REALNO  (VAR, 10, 'ABUT* , ‘ABUT’ ) 

55  CALL  GRAF  (0.  ,ANS,AN,-.5,.25,.5) 

DO  60  1=1, N 

IF  (IB.LE.  2)  X(I)=. 1*REAL(I) 

IF  (IB.GT.  2)  X(I)=.  3*REAL(I) 

60  CONTINUE 

65  CALL  CURVE  (X,Y,N,0) 

CALL  ENDPL(O) 

RETURN 

END 
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SUBROUTINE  PLOTI2  (Y,X.N,IFIG,TITLE$,TIT$,IB,IT,IFL,TL,SN,VAR) 
DIMENSION  Y(N),X(N) 

CHARACTER*80  TITLE$,TIT$ 

REAL  X,Y,STD 
INTEGER  IFIG 

999  WRITEC*,*)’ SELECT  DISPLAY' 

WRITE(*,*)'  1  =  FIXED' 

WRITEC*,*)’  2  =  ADJUST' 

WRITEC*,*)’  3  =  NORMALIZED’ 

READC*,*)IAD 
YMAX=1.  0 

IF  CIAD.EQ. 1)  THEN 
WRITEC*,*) 'SELECT  YMAX  CREAL)’ 

READC*,*)YMAX 
YMIN=-YMAX 
DO  10  1=1, N 

IF  CYCD.GT.YMAX)  YCI)=YMAX 
IF  CYCD.LT.YMIN)  YCI)=YMIN 
10  CONTINUE 
END  IF 

IF  CIAD.EQ.  1)  GO  TO  111 
YN=0. 

DO  15  1=1, N 

IF  CYCI).GT.YN)  YN=YCI) 

15  CONTINUE 

IF  CIAD.EQ.  2)  YMAX=YN 
IF  CIAD.EQ.  3)  THEN 
DO  20  1=1, N 
YCI)=YCI)/YN 
20  CONTINUE 
END  IF 

111  IF  CIB.LE.  2)  A^^=0.  1*REALCN) 

IF  CIB.  GT.  2)  AN=0.  3*REALCN) 

ANS=AN/4. 

C  CALL  PRTPLTC76,2) 

C  CALL  TEK618 

CALL  PAGE  Cll.  ,8.  5) 

CALL  NOBRDR 
CALL  HWROTC' MOVIE') 

CALL  AREA2D  C5. ,2. 5) 

CALL  HEIGHT  CO. 15) 

CALL  YNAME  C ' AMPLITUDE? ', 100) 

CALL  XNAME  C 'AZIMUTH  CMRAD)?' ,100) 

CALL  MESSAG  C 'FIGURE  $' ,100,1. ,5. 6) 

CALL  INTNO  (IFIG, 'ABUT' ,’ ABUT' ) 

CALL  MESSAG  (TITLE?, 100, 'ABUT' , 'ABUT* ) 

CALL  INTNO  (IB, 'ABUT' , 'ABUT' ) 

CALL  MESSAG  ('TARGET  NUMBER  $',100,2.48,5.4) 

CALL  INTNO  ( IT, 'ABUT' , 'ABUT') 

CALL  MESSAG  (TIT?, 100.2. 48,5.  2) 

CALL  INTNO  (IFL,  *ABUT% 'ABUT' ) 

CALL  MESSAG  ('TARGET  LOCATION  »  $’,100,2.48,5.0) 

CALL  REALNO  (TL, 1, 'ABUT' , 'ABUT' ) 

CALL  MESSAG  ('  MRAD  $' ,100, 'ABUT' , 'ABUT' ) 

CALL  MESSAG  ('S/N  =  $’,100,2.48,4.8) 

CALL  REALNO  ( SN , 4 , ' ABUT' , ' ABUT' ) 
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CALL  MESSAG  ('VARIANCE  =  S' ,100,2. 48,4. 6) 
CALL  REALNO  (VAR, 10, ' ABUT* , ^ABUT' ) 

55  CALL  GRAF  (0.  ,ANS, AN, -YMAX, 'SCALE* ,YMAX) 
DO  60  1=1, N 

IF  (IB.LE.2)  X(I)=.  1*REAL(I) 

IF  (IB.GT.  2)  X(I)=.  3*REAL(I) 

60  CONTINUE 

65  CALL  CURVE  (X,Y,N,0) 

CALL  ENDPL(O) 

WRITE(*,*) 'ENTER  1  FOR  ANOTHER  PLOT' 

READ(*,*)IPLT 

IF  (IPLT.EQ.  1)  GO  TO  999 

RETURN 

END 
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